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Calcium Hydroxide as a Highly Alkaline pH Standard 
Bower, and Edgar R. Smith 


Roger G. Bates, Vincent E. 
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1. Introduction 


The National Bureau of Standards standard pu 
e is defined in terms of several fixed points in 
guch the same manner as the International Temper- 
wre Scale. The primary standards are solutions 
hose pli values are only slightly affected bv dilution 
oe by accidental contamination of the solution with 
mees of acid or alkali from the walls of the container 
The substances from which 
de standards are prepared are, in turn, stable ma- 
wrals: thev are obtainable in the form of certified 
amples from the Bureau. The five standards thus 
established cover the pH range 1.68 Y.1S at 
mC! it is the purpose of this paper to describe 
be establishment of a sixth standard, a solution of 
alum hvdroxide saturated at 25° C, which will 
send the standard seale to pH 12.45 at 25° C 


gfrom the atmosphere 


to 


2. Need for an Alkaline pH Standard 


In practice, most pH values derived 
drectly or indirectly from the emf of the cell 


ure either 


solution X 
KC], He.Cl 


3.5 \J or 
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Pt; H, 4 or glass electrode 


sata 


where the vertical line marks a boundary between 
two liquid phases The transfer of positive and hega- 
live ions at different rates across this liquid junction 
gves rise to a potential difference. Inasmuch as the 
imnsference numbers of the ions vary through the 
boundary, as do the concentration gradients, the net 
targe transferred is rarely zero. Furthermore, the 
sgn of the potential difference may be either positive 
rhegative, and the magnitude can neither be meas- 
ured exactly nor calculated If included in’ the 
measured emf, this potential would result in errone- 
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Fortunately, the concentrated solution of po- 
tassium chloride that composes the salt bridge re 
duces the liquid-junction potential to a small con- 
stant value for most solutions of intermediate 
acidity, pH 3 to 11 When solution X contains iup- 
preciable amounts of the highly mobile hydrogen 
and hvdroxvl ions (pH less than 3 or greater than 11 
however, the liquid-junction potentials may differ 
considerably from the relatively constant values 
maintained in the region of intermediate pH. For 
accurate pH measurements in these regions of high 
acidity or high alkalinity, therefore, it is particularly 
important that reference standards of low and high 
pu be available. 

This need is readily demonstrated. A pu measure- 
ment is essentially a determination of the difference 
between pHx, the pH of the unknown solution, and 
pugs, that of the standard. If /y is the emf of cell 1 
and /’s that of the same cell when the standard is 
substituted for solution X, and /yy and /yg are the 
corresponding liquid-junction potentials, we have 


Ey.—Esz E EB 
~26RT F* 2.3026R7 


x 


He =. 
/ kf 


pus 9 
where R, 7, and F are the gas constant, the absolute 
temperature, and the faraday, respectively 

There is no wav to evaluate the last term of eq 
Henee, for the usual pH measurement the experi- 
mental conditions are chosen so as to nullify the 
greater part of the junction potential The two 
potentials, Aig and yy, are then assumed to 
nearly equal and to cancel rather completely when 
the difference is taken (eq (2 Hence the difference 
of pH. is assumed to be proportional to the difference 
of emf, Ey-Fs. Evidently the residual liquid- 
junction potential will be small when the standard 
solution and the unknown solution differ but little 
in pH, for the concentrations of the free hydrogen or 
hydroxyl ions (on which the liquid-junction potential 
largely depends) will be nearly the same in solution 
X as in solution 3, 


) 
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3. Calcium Hydroxide as a pH Standard 


The ideal pu standard is a stable solution easily 
reproducibly prepared from pure materials. 
The pH of the solution should not be markedly 
affected by changes of temperature. The 
value [6] should be high and the dilution value [7] 
low; in other words, the pH of the solution should 
not be altered appreciably by contamination with 
traces of acidic or basic impurities or by a change in 
the total concentration of the buffer substances 


The ionization of water is strongly affected by 
temperature changes, and hence the pH of most 
highly alkaline aqueous solutions is also sensitive to 


temperature. The pH falls) with 
by as much as 0.035 unit deg 
(" {8}. This imposes, of unfortunate 
limitation on the use of an alkaline pH standard 
It should be noted, however, that a certain degree 
of temperature control is necessary for pH measure- 
ments in the highly alkaline region, not only because 
the pH of the standard is sensitive to temperature 
“unknowns,”’ in general, 


alteration of 
increasing temperature 
course, an 


changes but also because the 
are similarly affected. 

The buffer value, or buffer capacity, 
containing appreciable concentrations of hydroxide 
ion is uniformly larger than that of buffer solutions 
of comparable concentrations in the pH range 3 to 
11, but these same solutions usually undergo a larger 
pH change on dilution [Ss] A high buffer Capacity 
is particularly important in an alkaline pH standard, 
for contamination with atmospheric carbon dioxide 
during and can never be completely 
avoided 

In principle, a solution of a moderately 
strong base will serve as an alkaline standard. In 
view of the high buffer capacity of aqueous solutions 
of even incompletely ionized bases at high pH, the 
addition of a salt of the is unnecessary. How- 
ever, it is very difficult, to find a 
strong base, either inorganic or that is a 


of solutions 


storage use 


strong or 


base 
if not im possib le. 
organic, 


stable solid The determination of the concentra- 
tion of the solution by a titration with standard 
acid is not only inconvenient and time-consuming 


but compounds the errors of three separate analytical 
Buffer solutions of trisodium phosphate 
same objections for, ow ing to the 
the phosphate 


operations 
are subject to the 
difficulty of preparing a sample of 
salt of the correct composition, they must be made 
from the secondary phosphate and standard alkali 
A pH greater than 11.2 at 25° C cannot readily be 
obtained with solutions of the alkali 
To be useful, an alkaline standard should have 
a pH of 12 or above. 

Substances that form highly alkaline solutions 
usually become contaminated gradually with car- 
bonate from exposure to the atmosphere when the 
container is opened. Many are also hygroscopic, 
and the moisture they acquire is not easily removed 
by simple drying procedures. If the carbonate were 
insoluble in water, however, the contaminant could 
be easily removed by filtration of the solution, and 
later contamination would be clearly evident when it 
occurred. If, in addition, the alkaline material were 


carbonates 


most 
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of moderate solubility, the saturated sol, 
serve as a standard. Under 
presence of carbonate would cause ho 
thermore, no weighings would be 
Because of its ease of preparation, 
solution of slaked lime, calcium hydroxid 
been employed as a reference solution of 
linitv [9, 10, 11). Tuddenham and Anderson fy 
saturated calcium chloride solutions with ; 
hvdroxide to achieve pH values from 11] 


On Wou! 


these con L1Ons, ¢] 


hecessary, 

Saturat, 
has oft, 
high alks 


Caleiy 
to Over 19 


as desired. The concentration of the Saturated gp. 
lution of calcium hydroxide in pure water js abo 
0.02 AM and the pH about 12.4 at 25 C A fox 
measurements made in an earlier study [12 howeve 


suggested that the saturated solution, like that of t} 
more soluble barium hydroxide, 


ciently reproducible to serve as a primary standap 


Later work has shown that the solubility of sample 


of calcium hydroxide free from soluble alkalies 
salts is indeed re producible within 1 percent, 

The solubility of calcium hydroxide decreases wy) 
rising temperature, being about 4 percent higher , 
20° than at 25° C, and 4 percent lower at 39 
Nevertheless, a separation of solid phase does po 
usually occur at 50° or even at 60° C when a Solutio: 
saturated at 25° is filtered and subsequently heate 
to those temperatures | 13, 14]. A solution of calein 
hvdroxide that has been saturated at or near 25° ( 
ean accordingly be employed asa pu standard ove 
considerable range of temperature 

The properties of a solution of calcium hydroxid 
saturated at 25° C are compared with those of ¢) 
phthalate and borax pH standards in the followin 
summary: 


. hvdr B 
, htha 
Molar concentration 0. 0203 0.05 0. 01 
Dens sitv (25° C 990] 1 OO17 
suffer value (3), mole 
pH unit ca. 09 0. O16 020 
Dilution value 
\pul , pH units 
for 1:1 dilution ca 2S 05 0] 
dp dt, pHunits degC 033 OO12 OO82 


It is seen that the pH change, in absolute measur 
on dilution of the calcium hydroxide standard is mo 
than 5 times that on dilution of the phthalate s 
lution of pH 4.02 whereas the temperature coeffi 
is 4 times that of the borax standard (pH 9.2 

appears that any highly alkaline standard will 
subject to these limitations in much the same degr 


4. Method of Assigning pH 


The emf method by which standard pH valu 
were assigned to solutions of calcium hydroxide | ; 
been described [1, 2, 8]. It will be summariz 
briefly here, but certain novel features of the preset 
treatment will be explained in detail 
lilution value attains large positive 
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The first step consisted in the measurement of the 
of of cclls with hydrogen and silver-silver-chloride 


trodes, containing the calcium hydroxide solu- 
ns W added potassium chloride: 
Pt: H. (¢), Ca(OH), (m,), KCI (m,), AgCl; Ag. (3 


moles per 
a solution does 


e svinbol m represents molality, or 
logran of water The molality of 
ot change With temperature and is therefore favored 
er volume concentration in measurements of this 
pe. several different molalities of calcium hydrox- 
le m ti range 0.015 to 0.0203 were selected. For 
sch calcium hydroxide concentration, three different 
wlalities of potassium chloride, namely, 0.015 m, 
01 m. and 0.005 m, were studied 

The second step was to compute values of the 
wantits pw, defined by eq }), for mixture 
f ealelum hydroxide and potassium chloride and to 
eermine, by extrapolation, the value of pwH 
e limit approached by pwH in calcium hydroxide 


each 


ylutions as the concentration of potassium chloride 


The value of pwH is computed 
and the standard 


ras reduced to zero 
rectly from the measured emf, £, 
wtential, /°°, of the cell [15] by 


KE 
2 3026RT FE 


pwH log -log Me, (4) 


luton 


shere f is an activity coefficient on the molal scale, 
is the gas constant (8.31439 } deg mole ,_ Fis 
e faraday (96493.1 coulombs equiv and 7 is 
etemperature on the Kelvin seale (deg C+-273.16) 
The values of pwll obtained by eq (4) for each of 

three concentrations of potassium chloride 
roved to be a linear function of log m,, the common 
garithm of the molality of calcium hydroxide, in 


e limited range studied (0.015< m,<7 0.0208 The 
mstants of the equation 
pw ah log my 5 


re determined by the method of least squares and 
sed to compute pwHl for five selected concentrations 

calcium) hvdroxide (m,=0.015, 0.0175,, 0.019, 
(2, and 0.0203) at each of the three chloride con- 
ntrations and at the 13 temperatures included in 
estudy. In agreement with our earlier observa- 
ons, pwil for a given value of m, was found to vary 
early with ms, the m ality ol potassium chloride. 
‘he limit, pw . in the absence of chloride was then 
eermined by fitting the values of pwH for the five 
vleeted caleium hvdroxide solutions to the equa- 


on 


pwHh = pwH 


Clits 


wv the method of least squares 
In the third step, the procedure necessarily departs 
irom thermodynamic rigor, for the pu lacks exact 


thermodynamic definition. The formal refitionship 


between pwH® and pHs is 

pug log fatty = pwH log fer, (7 

where fo, is the activity coefficient of chloride ion 

in the chloride-free solution of calcium hydroxide 

The NBS standard pH seale is based on the con- 

vention that the last term of eq (7) be expressed by 
the Debve-Hiickel equation 


A 
log I, a NX 
] Ba vi 


with a reasonable choice of the ion-size parameter, d 
In eq (8), wis the ionic strength, while A and FP are 
constants for the water medium at a particular 
temperature [16]. The most probable value of a, 
for chloride ion in a mixture of simple electrolytes 
is regarded to lie between 4 and 6, but values as 
low as 3 or as high as 8 cannot be termed unreason- 
able. Inasmuch as the convention adopted does 
not specify a single value of a;, the pHsg is not defined 
exactly through eq (7) and (8). The uncertainty 
amounts to a few thousandths of a pH unit at an 
ionic strength of 0.01 and to slightly more than 

0.01 unit at an ionic strength of 0.1 [8]. 

Because calcium hydroxide is not completely 
ionized,® the ionic strength of calcium hydroxide 
solutions is somewhat than However, 
the hydroxide ion molality, mou, can be evaluated 
with adequate accuracy from pwH® and the ion- 
product constant, A,,, for water at the appropriate 
temperature 


» 
OM). 


less 


log Mon log K t pwH {)) 
The first dissociation step is presumed to be complete 
and the second incomplete; hence, the molalities of 


the ions Ca** and CaOQH® are given by 


Wig Noy m 10) 
and 
Mean 2 mM, Mion, 11) 
and the ionie strength is 
M 2 Miou Mt). 12) 


5. Experimental Procedures and Results 


5.1. Solubility of Calcium Hydroxide 


There is disagreement among literature 
values for the solubility of calcium hydroxide, as 
well as some indication that the apparent solubility 
is dependent upon crystal size [18, 19]. Bassett [18] 


some 


See, for example, the paper of Bell and Prue [17 he second dissociation 
constant is about 0.05 
‘ Eq (9) depends for its validity upon the approximate equality of the activity 


ons of like charge in solutions of moderate concentra 
of water in such solutions does not depart 
wecurate the more dilute 


coefficients (fc; and fou) of 

mand upon the fact that the activity 
greatly from unity rhe formula is consequently more 
the solutior 
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wide an extensive study of the solubility from 
» 100° C and lists two values over a considerable 
of temperature; his ‘“‘true solubility’ is from 
to 15 percent lower than the solubility of very 
small crystals. The true solubility at 25° C, in 
oles per kilogram of water (molality), was found 
Bassett to be 0.02018, in reasonable agreement 
with 0.0205 found earlier by Moody and Leyson [20], 
but significantly lower than 0.0209 obtained by 
interpolation between the data at 10° and 42° as 
determined by Haslam, Calingaert, and Taylor [21] 
For the solubility of large crystals at 25° C, how- 
ever, Johnston and Grove [13] found a value of 
0.01976 m, and quite recently Peppler and Wells [19] 
have concluded that the solubility of large, well- 
formed crystals is 0.0184 m at 30 Inasmuch as 
the solubility of calcium hydroxide appears to be only 
about 4 percent lower at 30° than at 25°, there 
remains a discrepancy of about 2.5 percent between 
these two determinations of the solubility of large 
crystals 
In the present study, the solubility of 
sumples of calcium hydroxide prepared in various 
Ways determined by shaking the material 
vigorously with about 15 times its volume of water 
im a vlass-stoppered bottle Between periods ol 
shaking, the bottle Was maintained uta controlled 
temperature by immersion in a water thermostat 
The excess of solid was removed on a sintered-glass 
funnel of medium porosity, and the concentration 
of the filtrate was determined by titration with a 
standard solution of hydrochloric acid to the endpoint 
of phenol red. Weight burets were used in the titra- 
tions. The solubility determinations were repro- 
dueible to 0.5 percent or better 
In the initial experiments with caleium hydroxide 
of reagent crade, solubilities us high as O. OBPOU m alt 
25° C were obtained, but this figure was lowered 
somewhat by repeated digestion of the hvdroxide 
under hot water. This behavior could result from 
the extraction of soluble alkaline impurities or from 
a gradual stabilization of crystal form during the 
digestion. Thereafter, the samples of hydroxide 
were prepared from “CP” low-alkali calcium 
bonate and the effeet was no longer observed 
The finely granular calcium carbonate was ignited 
in platinum dishes at 1,000° C for about 45 min, 
and the resulting lime, after cooling in 
was added with stirring to water. The mixture was 
heated to boiling with continual stirring and was 
then filtered. The caleium hydroxide obtained in 
this wav Was dried at 110° C and powdered 
The solubility of five samples of hydroxide pre- 
pared in this manner was found to be 0.02037 m, 
with a standard deviation of O.00011, and this 
figure was substantially unchanged after seven ex- 
tractions of the product with water. To test the 
effect of aging, the caleium hydroxide was stored 
under water in a polvethvlene bottle and samples 
removed from time to time for solubility tests con- 
ducted as described After 1 month, the 
solubility of four samples wis found to be O O?032 m. 
with a standard deviation of 0.00006, and after 6 
months had elapsed a figure of 0.02022 m was ob- 


several 


Wiis 


i desiccator, 


above 


tained as the mean of two 
average solubility is 0.0203 mat 25°C. At 
solubility was found to be 0.0211 m and at 320 ( 
O.0O1L96 m. 

These figures are in excellent agreement with 4 
true solubility of Bassett {1S}, who found ( 021] . 
at 20°, 0.0202 m at 25 
appears that they 


represent satisfactorily 


caleium hydroxide, prepared from low-alkali caleiyy 
carbonate, with water at 20°, 25°, or 30° ( 


of the material. It should be noted that a chang 


of | pereent in the concentration of calcium hydroy. 


idle corresponds to about 0.004 in the pul of th 
nearly saturated solution. 


5.2. Electromotive-Force Measurements 


The caletum hydroxide solutions for the PT 
measurements were prepared by dilution Of near 
saturated stock solutions whose concentrations ha 
been established by titration with 
A weighed amount of potassium chloride was th 
added to each solution flask, or a portion of a 0.8 
solution of the salt introduced from a welch buret 
The dilution and the addition of chloride 1th an 
mosphere of carbon-dioxide-free nitrogen were 
complished with the aid of an arrangement yy 
similar to that deseribed by Bates and Acree 22 
The chloride was a bromide-free fus 
sumple prepared in the manner deseribed in 
earlier publication 235] 

The cells were rinsed, flushed, and filled in 
usual manner. Solution was admitted from 
solution flask, the cell emptied and flushed with pu 
hvdrogen, and the evele repeated before the fir 
portion of solution was admitted The preparaty 
of the electrodes has been described elsewhere ‘ 
Kach cell contained two pairs of electrodes, and 1 
duplicate Measurements were averaged 

Initial measurements were made at 25° C, thos 
from 0° to 30° C on the second day, and those fro 
20° to 60° © on the third day With these high 


alkaline solutions. y= | 


standard nT 


potassium 


the final measurements at 25 
when they were obtained, did not neree as clos 
with the initial values as those for solutions of | 
or intermediate pu often do Differences of 0.3 
0.5 mv were not uncommon. The emf data wi 
corrected in the usual way to a partial pressur 
1 atm of hydrogen 


6. Determination of pwH and pwH* 


A value of pwH was computed from each correct 
emf value by eq 1), and the constants @ and 5 of: 
(5) were obtained for each of the three molalities 
potassium chloride (m, equal to 0.015, 0.01, or 0.005 
in the cell solutions. The data are summarized 
table 1. The third column gives the number of e 
studied, and the sixth column lists o, the standa 
deviation of a single pwH value from the line defi 
by the constants da ana h The last five columns g 
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7. pH of Calcium Hydroxide Solutions of a, on the calculated pHs is shown in figure 9 9g, St 


The solid line passes through the values listed ir 


lhe ionie strengths of the five solutions of calcium | table 2, and the dotted Imes indicate the lisplace. FAs | 
ivdroxide at 0°, 25°, and 60° C, computed by eq | ment resulting from a choice of a4,;=3 and 8 ot of 
9) and (12), are as follows ‘raigh 
one = } 71 

Taspie 2. pHs of five solutions of calcium hydroxide frop, thet 

fo BO° ¢ itage 

wat ted 

2 ; { , te 

er © 25° © 60° C g to 
( prop 
1s l we 0 l r } ol 

0. 0205 0. O54 0. 049 0. O47 1 13 004 —-- 12. O7S = + 9 aaa ting 

0. O53 O49 O46 , 12 S00 12 Sf 12 THe 12 754 p~ il 7 
Ol O50 O47 O44 J 12. 620 12. 623 12. M02 12 ‘ 2 5K emp 

O75 O47 O45 O41 9 12. 454 2 449 12 429 12 a8 isiaa The 

OLS O40 O37 O85 4) 12 2M 12 20] 12. 27 12.2 2 174 the 
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The values of pli viven in table 2 and plotted 0) 11.704 11.701 11. 6S 11 ¢ 1 506 + no 
as a function of the molality of calcium hydroxide = +4, oe 1. on +t “4 rbon 
in figure 1 were calculated by eq (7) and (8) with the mat 
use of these ionic strengths. Each pHg value is the Saturated at 25° ( west 
mean between that obtained with a,;—4 and a,—6; vy al 
these two different estimates of pH, differed by 3 eas —_ H 
0.009 unit for the solution saturated at 25° (0.0203 or t 
m) and by 0.007 for the 0.015-m solution. In the nel 
light of the statistical probable error in pwH?® and ib ie 
the partially arbitrary selection of a,, the pHg is P eth 
assigned an uncertainty of 0.01 unit The effect , na o} 
a a be 
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§ Standardization in the Range pH9tol2.5 


As has been pointed out earlier in this paper, a 
ot of the emf of cell 1 as a function of pH is a 
waight line of slope 2.302627 F when the pH. is 
ther ”) low nor too high At high alkalinities a 
jtage Geparture or change in slope is to be ex- 
ted There will be an accompanying aberration 


the p retical pu scale near its upper end, operat- 


g to eld pu values that are too low It Is 
propriate to consider next the probable magni- 
le of the error and means of minimizing or elimi- 
ting by proper standardization of the pH 
sembl\ 
fhe error has its origin in an experimental defect 
the pH method, namely the changing liquid- 
nevlol potential with alteration of the coneentra- 
1 of highly mobile hydroxide ions. It is often 
+ noticeable below pu 11. but earlier work with 
rhonate buffer solutions |12]) has indicated that 


may be detectable at pu 10) This earlier study 
gested further that the magnitude of the discrep 
yat 25° Cis about 0.01 unit at pH 10, 0.02 unit 

and about 0.05 unit at pu 12.6 The 


pH 11.1 


or thus appears to increase in a fairly recular 
mner with the pH of the test solution 

lt is of some portance theretore, to cdlete rmine 
ether this apparent regularity extends to solu- 


The difference of poten- 
between hvdrogen electrode dipping in a 
ition saturated with calcium hydroxide at 25° C 
d another dipping in the standard 0.025 MZ phos- 
ate butter pls ) Was therefore 
termined \ cell two hydrogen- 
compartments [12] was used. Contact be- 
solutions Wis established through au satu- 
The observed 


ns of calcium hydroxide 


one 


Sth) ni 45) + 


Consisting ol 


troa 
en thre 
ted solution of potassium chloride 
28S \ corre- 


ce of potential wis 0 at 25 


onding to a pul of 12.42, for the solution of calcium 


roxide This solution was assigned a pHs ol 
1454 from measurements of cells without a liquid 
Lior compare table 2 The discrepancy ol 
2! pil 12.45 1s reasor ably consistent with 

ear progression of the error trom pil 9.18 the 


les in the uppel 
the prea tical pu 


error of about O05 


H of the borax standard), which 
error-free reciol ol 
) ‘ a ) wet , aT 

lo p 1? SS. where a 
24 


COTS ULETICE il 


ol Tre 


ims been found 


As a 


H OIS will 


standardized at 
reading betwee! 


pul ete 
vield i shehtly low 
und 45 


LZ standardized at pul 
vield i 


the 
will be 
hay t’ 


nnd a meter 
shehtly high reading in 
evidently 


Hence 


HO IS 


245 will sume 
errol 


1] rit 
jp 3 itil 


We 


1.0105(LY— Es ; 
2 3026RT F 


ply plig 17.08(E, | os at 25° C 14 


ich as the liquid-jur etion potential Is SOTHDe- 


the 


I usil 


it dependent upor structure of the liquid- 


liquid boundary, it is advisable to determine the 
error for the particular assembly that is being used 
Freshly prepared borax and calcium hydroxide solu- 
tions should be used for this purpose. If the error 
differs much from 0.03 or 0.04 unit, the numerical 
coefficients of eq (13) and (14) should be adjusted 
accordingly. The difference as indicated by a pu 
meter may be influenced somewhat by errors i 
scale length and in the temperature compensator, if 
these have not been calibrated. 

Because the pH meter reads pH units directly, it 
mav be convenient to apply a correction to the indi- 
cated pu, rather than to compute the correct figure 
by eq (15) or (14 from a measurement of the emf 
The temperature compensator of the meter is, how- 
ever, designed permit the selection of different 
values of the slope pH emf, corresponding to the 
value of F)(2.3026R7) at different temperatures 
Hence, this device permits pu to be read directly 
from the seale of the instrument, even though the 
functional relationship between the emf of the cell 
and the pH value varies rather widely, as it may do 
when the temperature of the cell changes 

When the temperature compensator indicates 25° C, 
the instrument converts emf differences into pH dif- 
the relationship ApH = 16.90AF; 


to 


ferences according to 


at 20 (’, the conversion is made according to 
ApH =17.19AF It is evidently possible, then, to 
compensate the liquid-junction error at the high 


end of the seale, permitting correct values at 25° ( 


to be read directly from the meter, if the tempera- 


ture compensator is set at 22° C, for at this temper- 


ature fF) (2.3026R7) is not far from 17.08 (compare 
eq (14 The values of this quantity at other 
temperatures are 

f 

10 17. SO 

lS 17. 49 

() 17. 19 

2. 16. OO 

0) 16. 63 


The alkaline error ol the class electrode made from 
Corning 015 glass is very much lower in solutions of 
calcium salts than in solutions of lithium or sodium 
a saturated solution of cal- 
cium hydroxide, this source of error may be of con- 
cern at temperatures from 40° to 60° CC. Many of 
the newer commercial electrode glasses contain ho 
calcium, and hence the voltage departures in calcium 
solutions are probably negligible 


salts 25 However, il 


9. A Standard Solution of Calcium Hy- 


droxide 


\ solution of calcium hydroxide saturated at room 
temperature is recommended as a pH standard for 
the highly alkaline range. A considerable excess of 
pure, finely granular hydroxide is shaken vigorously 
in a stoppered bottle with water at room tempera- 
ture The CTOSS CXCESS of solid is allowed to settle 
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niperature recorded to the nearest degree Cel- 
and the suspended material removed by filtra- 
with on a sintere l-glass funnel of 

jum porosity 
It has been recommended that a slurry containing 
xcess calcium hydroxide be used for standardization 
11 However, large errors are sometimes 
meurred in pH measurements of suspensions [26], 
and it seems best to rule out this uncertainty by re- 


suction 


pul DOsSCS 


moval of the solid phase 

For the saturation, a polvethyvlene bottle is very 
satisfactory; in fact, it may be found convenient to 
keep the solid continuously under water in a well- 
stoppered polyethylene bottle ready for final satura- 
tion and filtration when a fresh solution is needed 
Contamination of the solution with atmospheric car- 
bon dioxide prior to filtration is obviously of little 
concern, Contamination of the filtered standard 
solution renders it turbid and is a cause for replace- 
ment 

The temperature coefficient of solubility is nega- 
tive, but the solution supersaturates readily and no 
precipitation of solid is ordinarily observed at 60° C 
Unfortunately, however, the change of solubility with 
temperature is sufficiently large to require that the 
saturation temperature be noted and the pus values 
column 2, table 2) adjusted by the appropriate 
amount. These values for saturation temperatures 


of 20°, 25°, and 30° C are 
Ne she Hsia ‘ Corres ) 
tio Solubil 
mp ‘ oe “— - 
rature 20) 25 10 = 
( 
20 0. 0211 12. 64 12. 47 iz. 3] O. O15 
25 0203 12. 63 12. 45 12. 30 0) 
0) ~O195 12. 61 12. 44 12. 28 Olb6 


The calcium hydroxide should be prepared from 
well-washed calcium carbonate of low -alkali grade 
The carbonate is heated slowly to 1,000° C and 
ignited for at least 45 min at that temperature 
\fter cooling, the calcium oxide is added slewly to 
water with stirring and the heated 
boiling. cooled, and filtered on a sintered-glass funnel 
of medium porosity. The solid is dried in an oven 
and crushed to a uniform finely granular state for use 

It is advisable to determine the concentration of a 
saturated solution of one portion of calcium hydroxide 
prepared from each particular lot of calcium car- 
bonate, deseribed earlier in this papel The 


Suspension to 


iis 





(MM) 


the 
solution are practically identical, differing by » 


molality (m) and molarity iturat, 


0.3 percent. If the concentration of the solyti 
saturated at 2! 
presence of soluble alkalies is indicated. If calejy 
carbonate of a higher grade cannot readily 

obtained, the impurity in the available carbonat, 
hvdroxide should be extracted by careful Washi 


with water 
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Heats of Formation of Hexacalcium Dialumino 
Ferrite and Dicalcium Ferrite 


E. S. Newman and Roald Hoffman 


\leasurements were made of the heat of solution of 6CaQO-2ALO>-FeO, and 2CaQ0-. FeO 
HC]. 26.61LH.O (2.000 normal at 25° C From these data. together with ippropriate 
eat-of-dilution measurements and data taken from the literat ire the heats of formation of 


the two compounds were calculated 


1. Introduction 


The major part of the iron in portland cement clinker is found in combination with lime 
and alumina in the form of solid solutions between the compounds 6CaQO-2Al1,0 3-Fe,O, and 
2CaO-Fe,O,[1]... One such solution, brownmillecite, had been recognized earlier as an important 
constituent of portland cement [2]. Thorvaldsen, Edwards, and Bailey [8] determined the heat 
of formation of this material, and Lerch and Bogue [4] measured its heat of hydration. A 
general investigation of the thermochemistry of the compounds occurring in hydraulic cements 
is being ¢ onducted at the National Bureau of Standards The heats of solution of solid solutions 
of 6CaO-2AL,0,-Fe,O3 and 2CaQO-Fe,O 3, measured in a mixture of nitric and hydrofluoric acids, 
were given in & previous papet DI The present paper gives the heats of solution of these two 


compounds in hydrochloric acid and the heats of formation calculated therefrom. 


2. Experimental Procedure and Results 


2.1. Thermochemical Equations 


Heats ot formation are usually obtained by the addition ofa considerable number of ther- 


mochemical equations The equations must be written with the proper coeflicients to ueTeE 
With the storchiometry Of the reactions they describe and to fit the experimental conditions 
chosen After a few experimental conditions have been arbitrarily selected, others in subsequent 
experiments will be fixed by the relationships among the equations To find these conditions 

entire series of thermochemical equations must be written in detail For this reason, the 


ermochemical equations for the determination of the heats of formation of 6CaQO-2AlL,0,-Fe,O 
ane > ‘aQ I (). are riven here Kor completeness, the heats of reaction are also given, although 


some of the values were obtained as the experimental results of the work described hereafter 


The arbitrarily s ted experimental conditions were based on the volume of the heat-of- 


solution calorimeter and the solub lities of the s ibstances involved The conditions chosen 
were such that a sample weight of 1 g¢ of powdered sample would be dissolved in 600.0 ¢ of 
2000 N HCl] (HEL. 26-61,H.LO at 25° C Upon this basis the following thermochemical equa- 
tions were writte! Che reasons for their use will be apparent or will be discussed later. Unless 


otherwise noted, the equations are written [or isothermal reactions at 2S ( 


a. Heat of Formation of 6CaO.2AI.,0,;.Fe,.O 


The heat of formation of 6CaO.2ALO FeO from the oxides may be obtained by SumMmMIng 


thre equations 


M4.1HCl, 6CaCl fAICL,, 2FeCl, 216S80H,O)(soln)—-6CaQO.2AL0,.Fe,O,(¢) — 814.1 HCL, 
26.61H,O) (soln), A//,;=—402.41  keal 


HC aQt IP(HC], 26.61H.O) (soln) —-6(CaClL, 53.23H.O) (soln) +6H,OU All, 280.58 


keal ~- 
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6H.O() +6(CaCl,, 53.23H,O)(soln)—6(CaCl,-54.23H,O) (soln), A//, 0.06 keal 


Fe.0O.(e, 20° C)+6(HCI1, 200H,O) (soln, 20° C)—2(FeCl,, 601.5H,O)(soln, 20° C), ATT, 
24.74 keal. 

Fe.0,(e, 25° C)—Fe.0,(e, 20° C), AH, 0.124 keal. 
b( HCI, 200H,O) (soln, 25° C)--6( HCI, 200H,O) (soln, 20°C), Al7, 107.13 keal 


2(FeCl,, 601.5H,O) (soln, 20° C)—>2(FeCl,, 601.5H,O)(soln, 25° C) All;=107.79 keal 


2(FeCl,, 601.5H.O) (soln)—2(FeCl,, 81.34H,O) (soln LO40H.O, As/,=4.028 keal. S 
6( HCI], 26.61H,O) (soln) + 1040H,O-—-6(HCI, 200H,O) (soln), AZ7, 2.568 keal. y 
PALO. (e,a 12( HCI, 26.61H.O) (soln)—4(CAICI,, 81.34H,O) (soln), AZ io 119.44 keal. 10 


6(CaClL, 54.23H.O) (soln) —784.1( HCI, 26.61H.O) (soln)—(784.1HCI, 6CaClL, 21190H.O) (soln), 
All 0.29 keal | 


tH AIC], 81.34H,O) (soln 7T84.1HCI, 6CaCl, 21190H.O)(soln)—(784.1HCI, 6CaCL, 4A1CI 
21520H.O) (soln), As/,.=0.48 keal |2 

2(FeCl,, 8$1.34H,O 7TS4.1HCI], 6CaClL, 4AICL 21520H,O) (soln)—(784.1HCI, 6CaClL, 4AICI, 
2FeCl,, 21680H.O) (soln), As/ 1.35 keal 13 


?) 


The sum of eq (2) and (3) is 


HCaQ le I2(HCI], 26.61H.LO) (soln 0) CaCl. 4 23HLO soln), Al/;, PS0.64 keal 14 


The sum of eq (4) to (9), inclusive, ts 


FeOnle 6(HCI], 26.61H.O) (soln 2(FeCl,, 8$1.34H.O) (soln), A// 22? 74 keal 15 


These two equations (14 and 15) represent the solution of two of the 


oxide constituents of 
6CaQ.2AL0O,-Fe.0, in HCI], 26.61H.O 


The sum of eq (1) to (13), inclusive, represents the formation of 6CaQO.2ALO,-Fe,O 
from the oxides, 


| ‘aQ) ¢ 2ALO (a Keo » (\¢ ‘al PALO keod op & A// IS 2Y 


keal, lt) 
b. Heat of Formation of 2CaO-Fe,O 


In a similar manner, the heat of formation of 2CaO-.Fe.O 


from the oxides was obtained 
from the equations 


306. 2HC], 2CaCl, 2FeCl,, 8419H.O)(soln)—2Ca0.-Fe.O,(¢ S1I6.2(HCL, 26.61H,0 


« soln , 

All 107.40 keal 17 

2CaQ(e LHC], 26.61 HO) (soln 2(CaCl, 54.23H.O) (soln), A//)< 93.55 keal IS 
FeoOs(e H6(HICI, 26.61 H,O) (soln 2(FeClL,, 81.34H,O) (soln), Al/ 22.74 keal lo 
2(CaCl,, 54.23H,0) (soln) 306.2 (HCL, 26.61H.O) (soln)-—+(306.2HCI], 2CaCl, S256H.O) (soln), 
A/1,,=0.11 keal 19 

2(FeCl,, 81.34H,O) (soln 306.2HCI, 2CaCl, 8256H,O) (soln)—>(306.2HCI1, 2CaCl, 2FeCh, 
S419H.O) (soln), Al7/. 1.34 keal 20 


The sum of the above five equations represents the formation of 2C'aO.Fe,O, from the 
oxides 
2C aO(e Fe.Oy (ec 2C'aO-Fe,O,(c), Alls, 7.44 keal (2 
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2.2. Preparation of the Samples 


The samples of dicalcium ferrite and hexacalcium dialumino ferrite used in this investiga- 
on were the same ones used in the previous work [5], and their preparation has been described 
because in 8 vears of storage the ignition losses had become about 0.2 percent, the samples were 
ntered overnight at 1,200° C and reground for 2 hr in a power-driven mortar just before the 
beginning of the present heat-of-solution measurements. The heat-of-dilution experiments 
ere performed with solutions of calcium chloride, aluminum chloride, and ferric chloride 
prepared from analytical reagents. The manufacturers’ lot analyses were assumed correct, 
because the heat effects to be measured were small and minor errors in the concentrations of 
the solutions would have a negligible effect on the measured values 


2.3. Apparatus and Procedure 


a. Heats of Solution 


All experiments were performed in an air-conditioned room controlled at 25 +0.8° C and 
au maximum relative humidity of 50 percent. The heats of solution were determined using 
|-¢ samples in an isothermal-jacket calorimeter [6] containing 600.0 ¢ of HCl, 26.61,H,O. 
The temperature changes were measured with a 25.5-ohm platinum resistance thermometer, 
Mueller bridge, and a galvanometer with an 8-m optical lever. The sensitivity of the measuring 
svstem was such that a scale deflection of 1 mm corresponded to a temperature change of 
about O.000073 deg + The corrected temperature rises were calculated according to the 
method described by White [7, p. 42]. Time-temperature data were taken with the aid of 
a chronograph during the rapid-rise portion of the experimental period, and the area under 
the time-temperature curve was obtained by using the trapezoidal rule. The 1l-g samples 
dissolved in 30 to 40 min. The powdered samples were kept in the balance case, the tempera- 
ture of which was recorded, and appropriate corrections were made for the sensible heat intro- 
duced into the calorimeter with the samples, using 0.2 cal/g for their specific heat. The final 
calorimeter temperature wa8 always within a few hundredths of a degree of 25.00° C, and no 
further corrections were made to bring the measured tsothermal heats of solution to that 
reference temperature The thermochemical calorie of 4.1840 absolute joules was used. 
The calorimeter containing the solvent acid was calibrated electrically [14, p. 217] over approxi- 
mately the temperature range covered by the heat-of-solution measurements. The standard 
cell and standard resistors used in measuring the electric energy were calibrated by the Electric- 
itv and Eleetronics Division of the Bureau, and the input was timed using a chronograph 
and the standard seconds signals originating in the Mechanies Division. The resistance 
thermometer was used i the customary manner without converting resistance changes to 
degrees Celsius. Hence the energy equivalent was obtained in calories per ohm, and is approx- 


nately equal to calories per 9.901 deg C 


b. Heats of Dilution 


The heats of dilution were determined in a calorimeter consisting of a 1-pt wide-mouthed 
vacuum flask partially submerged in a constant-temperature water bath. The cork and 
upper end of the flask were protected from room-temperature variations by a double-walled 
copper cap dipping into the water, as suggested by White [{7, p. 152]. The liquid sample was 
held in a thin-walled glass bulb with a long stem. The bulb was placed beneath the surface 
of the calorimeter acid with the stem held ina clamp above the copper cap A thin glass rod 
was placed so as to rest in the bulb with the upper end protruding slightly above the stem, 


and a 40-cm glass tube of larger diameter was attached to the stem with a short piece of 


rubber tubing An S-¢ weight was held at the top of the tube. to which was connected a 
100-ml reservoir containing air at a gage pressure of about 0.3 lb/in’. To introduce the sam- 
ple, the weight was allowed to fall, and the Impact drove the lower end of the rod through the 
bulb. A small quantity of air passed through the tube until the pressure in the reservoir fell 
to about 0.1 Ib/in?. The pressure was then released to allow the stem of the bulb to fill to 
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This somewhat elaborate system was used to insure co 


liquid level in the calorimeter. 
The arrangement of the calorimeter required 


te and prompt introduction of the sample. 
iat the bulb be placed off center, and the stem frequently broke if the bulb was sim} 
pushed against the hemispherical bottom of the vacuum flask. On the other hand, the rod 
frequently punched only a small hole through the bottom of the bulb, and the air pressu 
was necessary to make certain that the sample passed into the calorimeter acid 
Three-hundred-seventy-gram quantities of HCI, 26.61,H,O or of distilled water were used 
The heat-of-dilution experiments were performed by using proportions of 


in the calorimeter 
the final concentrations should be identical with those 


acid or of water to solution such that 
in the heat-of-solution determinations The concentrations and quantities to be used In each 
experiment were calculated from the appropriate thermochemical equation listed in section 
a The calorimeter was calibrated electrically, using the same procedure and equipment 
The heat capacities of the solutions contained in the 


deseribed in the preceding section 
sample bulbs were estimated from the apparent molal heat capacities of the solutes, and the 
energy equivalent of the vacuum flask calorimeter and its contents was increased correspond- 
ingly Here again the heat effects to be measured were small, and small errors in the heat 
capacity did not affect the final results significantly. The final temperature of the calorim 
eter was always close to 25.00° C, and, as in the heat-of-solution measurements, no further 
corrections were made to obtain the heats of dilution at the reference temperature. Using the 


2.3 the bridge current was reversed 


platinum resistance thermometer described in section 2.3.a, 
each time the calorimeter temperature was recorded, This reversal, by doubling the observed 
deflection, doubled the sensitivity of the measurements. The very small corrected tempera- 
ture changes during the experimental period were calculated by graphical or analytical extra- 
polation of the time-temperature data observed during the initial and final rating periods [{14, 
p. 232 

The heat effect of breaking the sample bulbs was determined by performing heat-of-dilu- 
tion experiments in which the solution in the bulb was the same as that in the calorimeter 
Under these conditions, no dilution effects would occur The average value of this blank 
determination was subtracted from each heat-of-dilution measurement 


2.4. Experimental Results 


In table 1 are given the results of the calibration experiments for the heat-of-solution 
calorimeter, and of the heat-of-solution and heat-of-dilution measurements The measured 


heat effects of breaking a number of sample bulbs such as were used in the heat-of-dilution 
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( ‘periments are also viven, The observed heats of dilution were corrected for the average 
observed energy of breaking the sample bulb before being calculated to the basis of 1 ¢ of 


owdered sample as given in table | 


2.5. Heats of Formation 


The procedure for calculating the heats of formation of the two ferrites from the oxides 
has been outlined by the equations in section 2.1. The experiments and calculations repre- 


sented by the equations will now be discussed 


a. Heat of Formation of 6CaO.2AlI.O,-.Fe.O 


Equation (1) represents the experimental determination of the heat of solution of 
HCaO-2ALO,-Fe.O,, converted to a molar basis. The equation is written in reverse for addi- 
tion, the sign of the heat effeet being changed. The experimental results in table 1 give the 
value 402.41 + 0.58 keal/mole for the heat of solution The standard deviation of the average 
is given here and elsewhere in this paper 

Equations (2) and (3) together represent the measurement of the heat of solution of 6 
moles of calcium oxide in an equivalent amount of hvdrochloric acid The heat of reaction 
AJ], was calculated from heat-of-formation data taken from Cireular 500 of the National 
Bureau of Standards [13 Revised values [11| were used for the heat of formation of HCI, aq 
The data for CaCl, solutions [13] were plotted and AH, was calculated as the slope of the curve 
of the heat of formation versus 7, at 1 —54 moles of water per mole of CaCl, 

The heat of solution of FeO, in HCI, 26.61H,0 cannot be determined directly, because 
FeO, will not dissolve Thorvaldsen, Edwards, and Bailey 3| determined the heat of solution 
in HCI, 200 at 20° C by an indirect method. Equations (5) to (9) represent the recalculation 
of their result, given as eq (4), to the appropriate conditions for the present investigation, i. e 
the heat of solution at 25° C in HCI, 26.61H.O. Equations (5), (6), and (7) represent the 
heat effect of changing the temperature of the reaction to 25° C, The heat capacity value of 
24.7 cal deg-mole used in calculating A//; was taken from Kelley’s compilation of low-tem- 
perature heat capacities [9 Revised values [11] for the apparent heat capacity of HCI, aq 
were used in estimating the value 29.3 cal/deg-mole for the apparent heat capacity of HCI, 
200H.LO over the range 20° to 25° C. This value was used in calculating A//,. The apparent 
heat capacity of FeCl, 601.5H.O, used in calculating A//;, was determined by a plot of heat 
capacity versus the square root of the molality, calculated from the data of Kangro and 
Klugge {10 The value 51.4 eal/deg-mole at 20° C was obtained. It was estimated [11 
that the apparent heat capacity would become about 1 cal/deg-mole more positive for cach 
degree rise in temperature, and the average value —48.9 cal/deg-mole was used for the change 
from 20° to 25° C represented by eq (7) 

To obtain the value of A//,, the dilution experiment represented by eq (8) (written in 
reverse) was performed. Ferrie chloride solutions are subject to hydrolysis and the formation 
of complexes, and relatively few heat-of-dilution data are available [10]. A few drops of HC] 
were added to the ferric chloride solution to stabilize it The heat effect possibly caused by the 
presence of this small amount of acid was ignored, as Thorvaldsen and his coworkers [3] found 
that the heat of mixing of (FeCl, 600H.O+ HCI, 200H.O) with HCI], 200H,O was negligible 
The heat-of-dilution value ?.014+0.014 keal/mole was found as the average of five deter- 
minations. Equation (9) represents the dilution of HCI, 26.61 H,O to HCI, 200H,O. The 
heat effect of this reaction was calculated from the revised values of the heat of formation 


of HCl, aq. The sum of eq (4) to (9) gives 
Fe,O,(e) +6CHCL, 26.61H.,O) (soln)—2(FeCl,, 81.34H,O) (soln), A//,;—22-74 keal. (15 


Aluminum oxide, like ferric oxide, cannot be dissolved in HCI, 26.61 H,O. Its heat of solu- 


tion was calculated in an earlier paper [12]. The calculation was repeated by using the revised 


heats of formation of hydrochloric acid solutions, and A//;) was obtained 
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Equations (11), (12), and (13) represent a dilution experiment in which solutions of Cat 
AIC], and FeCl, were added in succession to a single calorimeter charge of HCI, 26.61H.0. 
The proportions and concentrations to be used were calculated from the equations, and the 


experiment was performed a number of times. The detailed results are given in table 1. The 
average values obtained were calculated to the basis of 1 mole of 6CaQO-.2 ALO,-Fe.O,, and are 
given as 0.29, 0.48, and 1.35 keal for A//,,, A//,, and Alf, respectively. The correspond 
0.017 keal, respectively. 


* 


standard deviations are 0.010, 0.017, and 
13), inelusive, gives eq (16), representing the formation o 


The summation of eq (1) to 
6CaQ.2ALO,.Fe.O, from the oxides, 


HCaO(c) +2ALO,(c, a Fe.O,(¢)--6CaO-2 ALO, Fe.O,(c), AH — 18.29 keal lt 


The standard deviation of the heat of formation, estimated from the sum of the variances o 


the heats of reaction entering into its calculation, is about 0.6 keal/mole 


b. Heat of Formation of 2CaO-Fe.O 


The heat of formation of dicalcium ferrite was calculated in a manner analogous to that 


HCaO-2ALO,-FeO,. Equation (14) represents the heat of solution of 2CaQO-Fe,O, written in 


The heat of solution, calculated from the data in table 1, was 107.40 — 0.068 keal mole 


reverse 
. represents the heat of solution of 2 moles of CaO 


Equation (18), which is one-third of eq (14 
in HCl, 26.61H.O, with a heat effect of 93.55 keal. Equation (15) has been discussed 


The last two equations, (19) and (20), represent the appropriate heat-of-dilution measure 


ments in which first CaCl, and then FeCl, solution was added to a single calorimeter acid 


The experiment was performed, in proportions of solution to acid calculated from 


charge 
the equations, and the values 0.11 and 1.34 keal per mole of 2CaO-Fe.O, obtained for A// 


and As foo, respectively The corresponding standard deviations are 0.008 and O.O11 keal, 


respectively. The summation of the equations indicated in section 2.1.b gives eq (21 


representing the formation of 2CaQO-F¢.O, from the oxides 


2CaOle FeO c ?2C aD FeO ¢), All, 7.44 keal 2] 


The standard deviation of the heat of formation, estimated in a similar manner to that of 


HC aQ) » ALO on ©) . is about 0] keal mole 


3. Summary 


\leasurements were made of the heat of solution of 6CaO-2ALO,-Fe.O, and 2CaO-Fe.O, in 
HC], 26.61HL,O (2.00 N at 25° C Heats of dilution of aqueous CaCl, AICI, and FeC! 


in the same acid were also determined From these measurements and from data in the 


literature, the heats of formation of the two compounds from the oxides were calculated to be 


18.3 and —7.4 keal mole, respectively 
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Response of a Sodium-lodide Scintillation Spectrometer to 


10- to 20-Million-Electron-Volt Electrons and X-Rays’ 
H. W. Koch and J. M. Wyckoff 


The 


response of 


a large-ervstal sodium-iodide 


spectrometer was studied for individua 


monoenergetic electrons extracted from a 50-million-electron-volt betatron operated betweet 


land 20 Mey The energy 


resolution with crystals 5 inches in diameter by 4 


inches in 


length and 5 inches in diameter by 9 inches in length reaches its optimum value for electrons 


of 6 Nie v. 


where 


the total width of the pulse-height distribution at half maximum is 4 percent 


Syntheses of pulse-height distributions for 11- and 19-Mev monoenergetic X-ray photons are 


made by using the eXpt rimental eleetron pulss 
-height 


compare favorably with the measured pulse 
Palma TuVves Svntheses provide one of the 
response funetions in this X-ray energy range 
1. Introduction 
The pulse-height distributions produced in total- 


sorption spectrometers by monoenergetic X-rays 
ow a narrow peak with small tail at low pulse 
ghits 1} These are therefore 
table for the accurate study of continuous dis- 
butions of X-ravs and for analyses of the shape of 
continuum However, an analysis of the spec- 
output in terms of the number of X-ray 
each photon energy requires detailed 
formation about the function of the 
ctrometer to monoenergetic X-rays or gamma 
rsin the energy range of interest 


spectrometers 


mete! 
ofons atl 


response 


For energies below 3 Mev, gamma rays from 
g-lived radioactive sources are available at 
ficiently small energy intervals to permit the 


irate evaluation of the spectrometer response. 
At energies above 3 Mev, the available nuclear 
amma ravs are few, weak in intensity, and generally 
suitable for the determination of response func- 
ms. For example, the relatively monoenergetic 
p.y)Be* gamma rays are a 17.6-Mev gamma ray 
ompanicd by a broad line at 14.8 Mey The 
ver energy line is sufficiently broad as to obscure 
the interpretation of the shape of the 
6-Mev pulse-height distribution. Similarly, the 
pyc and the Be*(any)C™ reactions provide 
j- and 4.43-Mev gamma rays, respectively, that 
weak in intensity, and experiments with them are 
sequently troubled by large relative backgrounds 
The response function can also be obtained by a 
mte Carlo calculation [2]. However, these caleu- 
are tedious and very time-consuming even 
en performed with the aid of automatic computers. 
difficulty arises from the multiple interactions 
itmust be followed in some detail at energies above 
Mev 
In order to provide some detailed data on the re- 
onse of a large-crvstal spectrometer for 10- to 20- 
ley monoenergetic X-rays, the response to mono- 
ergetic electrons was examined [3] as described in 
sreport, and then the pulse-height distributions for 


mewhat 


Hons 


s research was supported by the United States Air Force, through the 
Pt / 

[Scientific Research of the Air Research and Development Command 

res in brackets indicate the literature re nees at the end of this paper 


-height distributions 


The computed results 


distributions for 11.6- and 17.6-Mevy 


few procedures for predicting spectrometer! 

two X-ray energies were synthesized from these data 
Specific X-ray energy values of 11 and 19 Mev were 
chosen for the svnthesis on the basis of a planned 


matrix covering the 10- to 20-Mev region. These 
values are reasonably near experimentally 
measured monoenergetic gamma-ray peaks of 11.6 
and 17.6 Mev, and so the calculations lend them- 
selves to a comparison with the experimental results. 

The sodium-iodide ervstals in the semtillation 
spectrometer used in the electron part of the experi- 
ment were 5-in. diameter by 4 in. long, 5 in. diameter 
by 5 in. long, and 5 in. diameter by 9 in. long. 
Pulse-height distributions were obtained for electrons 
with energies between 1 and 20 Mev. The results 
gave the shapes of the distributions, the variation of 
the relative resolution (full-width at half maximum 
divided by the pulse height at the maximum) with 
electron energy, and the linearity of the pulse height 
of the peak of the pulse-height distribution versus 
electron energy 

The pulse-height distributions for monoenergetic 
X-rays were svnthesized for a 5-in.-diameter by 9-in.- 
long-crystal spectrometer from the electron results, 
This work, which is also described in this report, was 
motivated by ihe fact that cood predictions of pulse- 
height distributions produced by X-rays are difficult 
to obtain. By using the experimentally determined 
pulse-height distributions for electrons with energies 
above 6 Mev, one removes the necessity for following 
the electron through its detailed interactions, 
would ordinarily have to be done in a complete 
calculation 


also 


as 


2. Response of the Spectrometer to Electrons 


The electron beam was removed from a 50-Mev 
betatron by a pulsed electromagnetic extractor [4]. 
Individual electrons with fixed energies in the range 
1 to 20 Mev were directed against the center of the 
flat surface of the sodium-iodide crystal, as shown in 
figure 1. 

Electrons were available from the extractor during 
a time interval of approximately 0.2 uw see every 
5,500 psec. The electrons from the betatron should 
be monoenergetic to at least 0.1 percent over periods 
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of the order of minutes [5]. Energy straggling 


of tin 
of the electrons produces the largest uncertainty in 


she energy in the present experiment because of the 
three 0.0014-in. aluminum vacuum windows and the 
alum! im-Magnesium-oxide entrance window of the 
ervsti contaimer. The most probable energy loss is 
00 Key 5], and the energy) uncertainty is less 
than 100 kev 


Adjustments on a beam of 10° electrons/sec were 
made to direct the electrons down the collimator hole 
ind to set the focusing magnet current without the 
spectrometer, The electron-beam Intensity was then 
udjusted so that the counting rate with the spec- 
trometer Was limited to five electrons detected per 
second. Since 180 bursts of electrons come from the 
hetatron in | see, 
trons Curing any one betatron burst Wis limited to 
shout 3 percent of all eleetrons detected A picture 
of the electron-beam size at the crystal Position is 
shown in figure 2 
The individual light pulses in the sodium-iodide 
rystals were detected by the photomultiplier-am pli- 
described 1] The 


the coincidence detection of elec- 


fier arrangements previously 


mutput pulses from the linea amplifier were sent 
nto a gain-of-two amplifier and into a 30-channel 
lifferential analyzer Two separate runs at each 
electron energy were made in order to check, first, 


the over-all shape of the pulse-height distribution, 
ind. second, the detailed shape of the peak The 
two runs were separated and followed by calibration 
runs With a mercury-relay precision sliding pulser 
ind a radioactive gamma-ray source. The 4.43-Mey 
RabDia,Be the 1.12- 


were used as calibration 


gamma ravs from a source ol 
Mev gamma ravs from Zn 
sources 

The betatron energy control developed by Saunders 
6] in conjunction with the beam extractor arrange- 
nt was used to produce electrons with energies of 


2, 1.96, 4.04, 5.12, 6.18, 38.30, 12.2, 13.8, 15.3, 
7.15, and 17.9 Mev, respectively. The energy of 
the electrons entering the ¢« rvstal was assumed to be 
00 kev less than these values 

Typical pulse-height distributions for 1.7-, 
.O-, and 15-Mev electrons entering the 5-in ~liam- 
ure shown 


8.0-. 


ter by 5-in.-long sodium-iodide crystal 

gure o 

Results of the percentage resolution as a function 
of the incident electron energy are given in figure 4 
The resolution for electrons of energy less than 6 
\lev deviates markedly from that expected on the 
basis of the statistical fluctuations in the number of 
photons deteeted by the photomultipliers |] This 
vas due to the seattering and straggling in the erystal 
ntry window, which are important at low energies 
At about 6 or 7 Mev the resolution values are those 
one would expect from the statistical fluctuations. 
Because the svnthesis ol pulse-height distributions 
ior X\-ravs that follows depends On pulse height dis- 
tributions obtained for liberated on the 
erystal axis and within the crystal, the energy reso- 
lution values used for electron energies less than a) 
Mev were those assumed to be determined by sta- 
tistical fluctuations with a Gaussian half-w «lth in- 
ferred from previous work with a similar speetrom- 


elect rons 
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ene t ete ed DV gamma-ray energy calit 
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he " ete rhe ver il line ibove the experiment 
e the t t il ment of energy t compar 
ny e-height p e elect d butior iwnd the low-ener 
m 
eter 1] For electron energies above 6 Nev. the 


experimental distributions similar to those in figure 
3 were used. 

Figure 5 proy «les a crude test of the linearity of 
the spectrometer response to electron energy. The 
ordinate is the kinetic energy of the electrons The 
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NUMBERS 


RELATIVE 


<1 scale is the potentiometer setting in volts 
betatron energy control [6] to which the elec- 
nergy has been shown to be related by a linear 


on The experimental points are the energy 


es assigned to the pulse heights of the peaks in 








crystal. No significant departures in distr: butioy 
shapes or resolution values were found for the differ. 
ent sized crystals. It was found that only 10 percep 
of the 15.6-Mev electrons produced a contriby. 
tion in the last 4 in. of the 9-in.-long crystal arrange. 

































ire 4 plus 300 kev. This assignment was made | ment. This contribution consisted of very small 
, comparison with a known gamma-ray energy | pulse heights. It was inferred from these testis thay 
143 or 1.12 Mev the energy loss from the three sizes of erystal used 
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3, Synthesis of X-Ray Pulse-Height Distribu- 


tions 


3.1. Pair Electron Contributions 


The most important contribution to the pulse- 
wight distribution of X-rays above 10 Mev results 
he energy loss in the crystal produced by the 


from 
pair electrons, Specifically, at 11 Mev, 68 percent 
nd at 19 Mev, 84 percent of the initial interactions 


rill be by the pair production process 

The origin of the electron pairs produced in the 
sin. by Q-in. ervstal will be distributed along the 
xis of the crystal if the A-ravs are assumed aper- 
nto a narrow beam that is directed along the 
xis. If it is further assumed that the electrons and 
ositrons behave sufficiently identically in their seat- 
ering and energy-loss characteristics, except for the 
voduction of annihilation radiation by the positrons, 
he resulting pulse-height distributions can be pre- 
ited by folding together the measured distributions 
lescribed im the last section for those electron ener- 
could combine to form a pair. The com- 
would be the incoming X-ray 
nereyv thins twice the ele tron rest-mass energ 


ined eleetron energy 
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There are obvious differences between the test 
electrons from the betatron and the pair electrons 
produced by an absorption of an X-ray photon 
One difference is the angular orientation of the elec- 
trons with respect to the crystal axis. If either of the 
two electrons of a pair has an energy less than 6 Mev, 
the probability that it will lose energy by brems 
strahlung is small. Therefore, the knowledge of the 
angles into which the electrons are emitted is unim- 
portant, because the electrons will lose energy by 
collision loss and produce i Craussian-type pulse- 
height distribution fig. 4 and its 
above On the other hand, if the electron energies 
are greater than 6 Mev, the electrons will be projected 
into angles less than 10 \7| For the present purposes, 
these angles are small enough to allow one to assume 
production of all pair electrons parallel to the erystal 
aXis 

Another difference between test and pair electrons 
is the crystal depth to the point of initial interaction 
The test electrons interact as soon as thes enter the 
ervstal, whereas the pair electrons interact after the 
depth of X-ray absorption. However, for the 5-in.- 
diameter by 9-in.-long erystal, 85 percent of the 
photons interact in the first 5 in. of the crystal. Also 
very little effect of length from 4 to 9 in. was de- 
tectable in the electron experiment described above 
Therefore. there should be little difference in the 
actual available interaction distance for test and pai 


see discussion 


electrons 

With the above assumptions, the pulse-height 
distributions produced by pair electrons were 
tained by folding together the pulse-height distribu- 
tions for electrons whose combined total energy 
equaled the X-ray energy of interest minus | Mev 
The individual distributions were inferred from the 


ob- 


experimental data such as shown in figure 4, and 
are drawn as histograms in figure 6. The results of 
folding together distributions of electrons with 


kinetic energies of 1 and 9 Mev, 3 and 7 Mev, and 


5 and 5 Mev are given in figure 7, a. Similarly, 
distributions of electrons of 1 and 17, 3 and 15, 5 
and 13. 7 and 11. and 9 and 9 Mev are given in 


figure 7, b. As all pair electron distributions are 
equally probable [S], these curves are simply added 
together, with double weight being given to the 5 
plus 5 Mev and the 9 plus 9 Mev curves, 
tively The annihilation radiation contributions are 
not shown in these curves but are treated separately 
in part 3.4. These contributions to the final syn- 
theses of the X-ray pulse-height distributions will be 
shown in figures 9 and 10. 


respec- 


3.2. Compton Electron Distributions 


The mMcoming X-rays with energies between 10 
and 20 Mev can enter the crystal and undergo 
Compton seattering as well as produce pairs. In 


order to include these contributions in the X-ray 
pulse-height distribution, it has been assumed that 
the distribution of electron energies for the initial 
Compton interaction of high-energy photons repre- 
the distribution of the sum of the electron 
energies after multiple Compton interactions of the 


sents 


same initial photons 
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The justification for this assumption comes from 
a consideration of the \lonte Carlo predictions ol 
Berge and Doggett 2) One ol their results {)] Wis 
the predicted distribution of electron energies result- 
ing from multiple interactions of an X-ray photon 
The histograms in figures Ss, a ana &, bare then data 
before a Craussian distribution has been added, lor 
an X-rav energy of 2.6 and 4.4 Mev, respectively 
Also drawn on these Saimie figures ure the Compton 
electron distributions for a single interaction [10] for 
the same X-rav energies. The similarity between 
the single- and multiple-interaction distributions is 
striking. As the X-ray energy increases, these 
distributions will become more peaked and more 
similar 

Because the Monte Carlo caleulation and = the 
single-interaction distribution are similar, the theo- 
retical distribution [10] was used to represent the 
electron energy distribution resulting from the Comp- 
ton seattering of X-rays in the range from 10 to 20 
Mev. This should be a valid assumption especially 
because the (‘ompton process is a correction and 
does not predominate in the determination of the 
shape of the pulse-height distribution for X-ravs 


3.3. Compton Scattering Electron Contribution to the 
Photoline Energy 


Those Compton seattering events for which the 
scattered X-ray interacts several times within the 
erystal, can lead to total absorption of the primary 
\-rav energy. These events have been included as a 
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Gaussian-shaped contribution at the photoline en- 


ergy The photopeak itself, arising from photons 
that have lost all their energy bby the photoeleetri 
process, is negligibly small at these energies.) hh 


the 10- to 20-NMey energy range, the area of this 
Craussian was related to the Compton electron dis- 
tribution by the following extrapolation from lowe: 
energies 

The photofraction, . has been defined yA as the 
ratio of the area under the photopeak of the pulse- 
height distribution. which Is represented by 
Gaussian, to the area under the entire pulse-height 
distribution Two constants characteristic of 4 
particular reometry ana energy enter the eCXpressiol 
for this ratio. The constant @ is that fraction of 
the total pair production events for which bot! 
annihilation quanta are absorbed im the crystal 
The constant 4 is that fraction of the (Comptot 
scattered photons that are multiply scattered an 
totally absorbed in the ervstal In terms of the 
photoelectric (7), pair (@), Compton (¢), and total 
vi absorption CTOSS SeCTIONS, p* is defined as equal LO 
T awn he u 

In ordet to extrapolate values of h to an energy ol 
19 Mev, 6 was ealculated up to 5 Mev from th 
data for a ana p* from reference 2 These cealeula- 
tions are shown in table 1 for a 5-in. diameter by 
9-in.-long sodium-iodide crystal. As will be observed 
from the sixth column, 6 varies slowly and should 
continue to vary slowly up to 19 Mev for a 5-in. bj 
9-in crystal As the Compton cross section (a) Is 
small and thus the product bo is small, errors du 
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the extrapolation of 6 in the evaluation of the 
-Nev photofraction values should be small 

The extrapolated values of 6 are given in table 2 
wether with the modified photofraction, ?*, from 
Mev, and the absorption coefficient data 
responding to those in table | Fr? is differenti- 
the latter values 
compared 


to 19 


ted Irom the low emerges p* because 


above 5 Nev eannot be 


r energies 


rectly with experimental values. 2?* contains a 
mponent sia, contrasted to am of p*, which includes 


total under the electron contribution 


nd not just the pair electron area under a Gaussian 


ured peal 
irve at tlre photoline POsILLON 
In svnthesizing the A-rav pulse-height distribu- 
on. thre Compton senattering events thraat result inh 
tributions at the photoline energy were imeluded 
sa Crussian wit} dn area h anid il half-width ol 
00, / PP The 
I Wiis 


urea under the ( ‘ompton electron 


stribution likewise tuken as I h 


3.4. Contributions Due to the Escape of Annihilation 
Radiation 


influence on the broadening of 
X-ravs will result from the 
scape ot one or two antl thilation photons from the 
the electron pair originally produced by 
\-ray Jerger and Doggett [2! have 
aleulated the escape of annihilation photons Irom a 
-in.-diameter by 9-in.-long found that 

} percent of the photons will be completely cap- 
pe reent of the events have exad tly 


The only 


e distribution due to 


remaiwine 


ositron of 


® comin Y 
crystal and 


red, and 21.5 
me annihilation photor escape with no interaction 
The remaining 15.2 percent of the events have some 
interaction of the one or two annihilation photons 
vfore some of their energy escapes from the crystal 


ror the present purposes it Is sufficiently vuecurate to 


assume that 36.7 percent of the events have exactly 
one annihilation escape and, therefore, that percent- 
age of the folded electron-pair pulse-height distribu- 
tion will be displaced downward in energy by 0.4 


ol 
Mev 


3.5. Synthesis of the X-Ray Pulse-Height Distributions 


The four contributions discussed above have been 
used to synthesize the pulse-height distribution for 
l1- and 19-Mey X-ray photons The areas of these 
contributions were all referred to the area under the 
folded electron-pair distribution. The relative con- 
tributions are given in table 3 

The component distributions are given in figure 9,a 
for 11 Mev X-rays and in figure 10,a for 19 Mey 
X-rays. The sum of distributions and the 
final result of the predictions are given as the solid 
lines in the 6 curves of these figures 


these 
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4. Discussion 


The svnthesized pulse-height distributions for 11- 


Niles and 19-Nl ey \-ravs were compared with (X- 
perimental distributions obtamed previously with 
11.6- and 17.6-Mev gamma ravs |1 


Figure 9.c, shows the experimental distribution for 
the B'(p.y) gamma-rays [1 
of this source was small, the channel window siz 
the differential pulse-height analyzer had to be made 
The « figure Ob, smeared 
with thre channel window used it the 
order to arrive at the prediction shown as the solid 
irve of The experimental and predicted 

irves have been normalized at the peak The ex- 
perimental curve is higher than the predicted one, at 
heights smaller than the one at which the peak 


Because the intensity 


Or 


lara urve of has been 


ball 


r xperiment Ith 


hnvure Ye 


occurs, for practically all pulse heights. However, 
the differences can be attributed to the Compton 
scattered gamma ravs trom the thick copper target 
issemblv, and to the collimator wall effects used in 
the work deseribed in reference |] The wall of the 


umetet collimator hole subtended a relativels 


=}! -( 


Sold al ol for the camma ravs trom thie target 





Hence gamma rays could interact and lose energy jy 
the wall and then penetrate into the large-diamete 
crystal. Both effects of target and collimator wou), 
result in degraded-energy gamma rays entering th, 
ervstal. Similar comments can be made revardino 
the comparison in figure 10,b, of the 19-Mev predje. 
tion with that obtained with the 17.6-Mev gamm, 
ravs from the Li‘(p,y) reaction [1] 

No detailed calculation to correct the experimentg 
curves has been made because the experimental daty 
(1), particularly that of the 11.6-Mev line, were to, 
Warrant extensive interpretations Never 
the comparison between prediction and ey. 


poor to 
theless, 
periment ts considered to be vood 

The method of svithesizing pulse-height distriby 
tions due to X-rays, which has been outlined, is th, 
best procedure known to the present authors in thy 
X-ray range from 10 to 20 Mev. Because of thy 
general and the expected validity of thy 
predictions, it is planned to apply a similar procedup 
for energies from 5 to 19 Mev in 2-Mev nmerements 
These results, combined with the work of Berger an 
Dovgvett 2, 
mental pulse-height distributions in the range frop 
100) kev 20) Alev, 
trometer with aio 5-in 


SUCCESS 


to 
diameter and =a 9-in.-loy 


ervstal are used 


Fuller IT) obtan 
for the 


Ps tree and I< (; 
by nin 


The work ol B 
mg a well-controlled electron 
experiment ts vratefully acknowl dved 


presel 
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leer 
me Hydrogarnet Formation in the System 
x th 
rdine . . > 
edie Lime-Alumina-Silica-Water 
Mims 
, Elmer T. Carlson 
Chia 
data 
Pe toy Glasses of the composition of grossularite and similar glasses poorer in silica were treated 
hvdrothermally, and the composition of the resulting hydrogarnets was estimated from the 
eve init cell size as calculated from X-ray powder patterns Che silica content of the hydro 
eX- garnets increased continuously with the temperature of formation approximating that of 
grossularit it SOO® C When mixtures of lime, alumina, and siliea, either as such or ir 
riby. Various states of combination, were treated in the same manner as the glasses, equilibrium 
: Was not attamed even after several weeks. At times as manv as three distinet hvdrogarnet 
su compositions were found to coexist Phe equilibrium temperature for the hydrothermal 
1 th decomposition of 3CaO-ALO>6HLO to 4CaO-3 ALO; 3HLO plus Ca(OH), was found to be 
L othy between 220° and 226° C Hydrogarnets containing increasing amounts of silica underwent 
th, the same tvpe of decomposition at progressively higher temperatures 
du 
ents 1. Introduction Such conditions are likely to exist in actual practical 
* ai processes involving formation of hydrogarnet ; hence 
peri The relationship between the isometric tricalcium it study of these variations should be of some interest 
ron ium ute hexahvadrate nnd evrossularite vrurnet was Finally there is the question of the upper limit 
spec- winted out by Flint Nile Murdie, and Wells | 2] of stability of the hvdrogarnets under hydrothermal 
long lt Wiis shown that these two compounds are the end conditions hint and MieMurdie, ana Wells 1] 
nembers of a series of solid solutions. of garnet-like | stated that a sample of tricalcium aluminate hexa- 
rystalline structure which were termed “hydro- hvdrate contaming a small amount of silica was 
ta varnets.” The general formula may be written stable up to 300° C, at which point it decomposed 
CaQ- ALO. mSiO.-(6—2m)HLO. where m mav vary to ealetum hvdroxide and gibbsite. However, it 
ma fom 0 to 3. The present work is an extension of | Was later shown by Johnson and Thorvaldson [4] 
if these earlier studies, undertaken for the purpose of | that pure tricalcium aluminate hexahydrate undet 
” nswering certain questions hot previously covered these conditions breaks down to form a less basic 
ompletely thuminate, 4CaO-3ALO,-3HLO, plus caletum hydrox- 
Th first of thes questions to be con sidered is thre ile This Wiis confirmed by Peppler and Wells 5), 
elation between temperature ana Dressure and the who placed the decomposition temperature at about 
aie omposition of the hyvdrogarnet in equilibrium with | 215 + A study of the effect of silica on this reac- 
. gs the fluid) phase Flint and Wells prepared their tion appeared desirable 
rdrovurnets Tromp mixttires of Varro‘ raw materimils The work of Flint, \MIleMurdie, and Wells included 
99. he desired lime-alumina-siliea ratios. selecting | also a series of quinary compositions, the fifth com- 
' l-and-error methods the optimum temperature | ponent beme ferrie oxide. This is omitted from 
| for each preparation. In general, the results | consideration in the present paper 
‘ Dlaiied al other Leribperatures Were Tol reported 
has been pointed out by Yoder [3] that theoreti- 2. Materials, Apparatus, and Procedures 
| lly the composition of the solid-solution phase at 
¥ ulibrium aaah be a funetion of Leriperature anid 2.1. Materials 
— This was supported by his data, on the The alumina used was a commercial preparation 
asis of which he constructed a pressure-temper- | of gibbsite (ALOs3H.LO) containing about 0.30 per- 
ture diagram showing the variation in refractive | cent of NasO, othe impurities being negligible in 
dex of che hydrogarnet phase With the exception | smount All other raw materials, except as other- 
Fone point, all his data were obtained at 500° C. Wise noted, were reagent grade chemicals. Ordinary 
higher. It was considered desirable to obtain | distilled water was used without further purification 
sme data at lower temperatures 
Another objective of the present study is to de- 2.2. Apparatus 
rmine how the course of the reaction is affeeted by Most of the reactions were carried out in stainless- 
chore of raw materials In theory, the composi- steel pressure hombea of the Morey type The 
mam OF tox procme: = & pen eenperasure and eeiet bombs were heated inh vertical tube furnaces. and the 
sure should be independent of the previous state of temperature was automatically controlled and re- 
a ~_~ map nevertheless it was found by lint corded by conventional upparatus using 1ron-conh- 
ind Wells that different results were obtained with stantan thermocouples. The temperature gradient 
agers starting materials. It must be assumed in each furnace was carefully explored, and the bombs 
i this was due to nonequilibrium conditions were positioned so that a thermocouple in the normal 
recording position and another at the bottom of the 
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amber gave readings agreeing within | deg. 
> minimized an error that might otherwise have 
as great as 10 deg, judging from preliminary 


No provision was made for independent pressure 
ontrol, hence the pressure developed in the bomb 
vas a function of temperature and of the degree of 
filling of the bomb. For temperatures below the 
eritical point, excess liquid water Was present, and 
the pressure was estimated from steam tables. For 
higher temperatures the pressure-volume-tempera- 
ture data of Kennedy 16] were used 

A cold-seal hydrothermal apparatus of the tywpx 
designed by Tuttle [7] was also used, permitting ex- 
tension of the study to higher temperatures. In this 
apparatus, both pressure and temperature were con- 
trolled manually A pressure of 15,000 psi was used, 
except for a few experiments on the effect of varving 
pressure 


2.3. Procedures 


The solid reactants, lime, alumina, and silica, may 
be introduced in a variety of forms or states of com- 
The separate oxides, or hydrates, for ex- 
mixed. When this procedure was 


bination 
ample, may be 


used. the lime was usually present ns the oxide 
freshly prepared by calcining the carbonate), the 
alumina as gibbsite, and the silica as gel (“silicic 


acid” 

Alternatively, two or more of the components may 
be chemically combined. For example, tricalcium 
aluminate was prepared by heating a mixture of 
calcium carbonate and alumina; likewise, tricalcium 
silicate from caleium and Tri- 
calerum aluminate hexahvdrate Wiis prepared i 
by hydration of tricaleitum aluminate at 150° to 
200° C. and (b) by the reaction between lime and 
freshly precipitated alumina at boiling temperature 
The latter method was also modified by the 
cipitation of small amounis of silica with the alumina 
Three 


this method 


carbonate silica 


copre 


hvdrogarnet preparations were obtained by 


contamime O10. O26. and O.5S mole of 


S10. per mole of ALO respectively However. as 
is shown later. in the last of these not all the siliea 
was present im the hvdrogarnet In leu of a svn- 
thetic aluminum silicate, Georgia kaolin, activated 


bv heating at 750° C, was used for a few experiments 

Finally, there is the possibilitv of using ternary 
glasses in which the three components are present in 
the desired proportions Cilasses were prepared by 
heating mixtures of calcium carbonate, gibbsite. and 
& Globar-type furnace The alu- 
mina ratio was kept at 3 to 1, whereas the silica 
alumina ratio ranged from 0.4 to 3.0 \iixtures that 
did not fuse at the highest temperature of the fur- 
about 1,450° C) or that devitrified on cooling 
were given further treatment 
In an oxvgen-gas blast until fused to a clear glass, 
The 


before 


silica gel in lime 


nwace 


Fragments were held 


and the molten drops were quenched ih Water 
class vround to puss a No. 200 
hvdrothermal treatment 
The starting materials, 
gram, usually were mixed 
slurry in a small platinum crucible, which was then 


Wiis s1eve 


totaling a few tenths of a 


with water to a paste or 





placed in the bomb with the desired quancity of 
water and kept in a furnace at constant temperatyy 

+ 1} deg C) for several days. When the a ld-seg] 
apparatus was used, the weight of solids averaged 
about 0.05 @, and the amount of water was & 
between 0.02 and O.05 2. The mixture was ¢ Closed 
ina small platinum capsule, which was then placed 
in the bomb; the bomb was connected to the pressyp 
line, the pressure raised, and the heat applied. Th, 
warmup period averaged about 30 min. After g 
heating period ranging from a few hours to 3 days 
the bomb was cooled in air (usually requiring aboy 
15 min) and opened, 

As a general procedure, the solid reaction produets 
were filtered off and dried by washing with alcohol 
and ether. Frequently, when the Morey bombs wer 
opened, the reaction products were almost dry; in 
such cases the drying was completed in an oven at 
about 100° C, without intervening filtration 

The dried products were first examined under thy 
polarizing microscope, after which X-ray powder dif. 
fraction patterns usually were obtained. The pat. 
terns were made with a recording X-ray diffractom. 
eter, using copper Ka radiation The X-ray pat- 
terns served not only for identification of ervstallin 
phases, but also as a means of estimating the compo- 
sition of the hvdrogarnet solid solution It was 
shown by Flint, MeMurdie, and Wells [1] that th 
hydrogarnets crystallize in the isometric svstem, an 
that the unit cell size is a function of COMPposi tion 
In the solid solution series $CaQO-ALO,-3510, 3Ca0 
ALO,-6H.O, the cell constant found to rang 
from 11.84 to 12.56 A On the basis of data from 
Yoder [3] a diagram show- 
ing the relation between composition and cell con- 
stant im this The relation is shown to by 
nearly, but not quite, linear. Thus if the cell con- 
stant is known, the COMposilion can be estimat 
with fair precision 

The cell constant 
the interplanar spacings, 
from the X-ray. diffraction 
measurement of the angles it 
sharp, distinet pattern, also to have a standard pres- 

In many of the reaction prod- 
the following 
~ present, serving as 


erally 


Wiis 


several sources presents 


SCTLOCS 


readily calculated {ror 
turn are derive 


mav be 
Which in 
angles Kor accural 
Is necessary to have 
ent for comparison 


uets discussed in sections, calcein 
hvdroxide wa 
standard. The caleulated unit cell sizes are believe 
to be accurate to 0.01 A for the determinations | 


which an internal standard was present, and + 0.05 A 


il reliable interna 


otherwise 

In numerous cases 
clear ey idenee of the presence ol Iwo, or occasionally 
three, separate hvdrogarnet COMpositions ih the 
same reaction product. It that, 
such cases, equilibrium had not been attained. Al 
though the condition usually could not be detected 


as is discussed later, there was 


seems evident 


microscopically, it was disclosed by the X-ray pat- 
terns, in which many of the prominent lines appeared 
either as doublets or as pairs (occasionally triads) ol 


Wher 


lines, depending on the relative COMpositions 


this occurred, only the low- and medium-angle reflee- 


tions could be relied on. as there was considerable 


overlapping at the higher angles 
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3. Results and Discussion 


3,1. Relation Between Temperatuie and Composition 


The phase equilibrium relationships involved in 
this study are complex, because the system is qua- 
the temperature range extends bevond the 
point for water, and the phase of primary 
solution A 
theore tical aspects Is bevond the scope of this paper 
However, it may be mentioned that there can be as 
many 2s four solid phases (five, above the critical 
temp ature) existing together at equilibrium, hence 
the presence of other solid phases along with the 
wdrogarnet is not evidence of failure to reach equi- 
Qn the other hand, the coexistence of two 


tera \ 
criti i 
terest Is a 


solid discussion of the 


ibrmum 
members of the same _ solid 


lack of attainment of 


solid that are 
solution series would indicate 
quilibrium This observed ina 
of instances reported in the tables The 
that hvdrogarnet phase indicates 


phases 
condition wis 
number 


onverse a single 


Ni vert! eless 


equilibrium, is not necessarily valid 


in the absence of any better criterion, the observa 
tion of a single well-defined hyvdrogarnet phase 
independent of the duration of treatment within thi 
limits of the experimental conditions, is herein taket 


to indicate the attainment of equilibrium 

After numerous experiments with various starting 
materials it was concluded that the glasses vielded 
equilibrium products more readily than did any of 
the other combinations. The data in table 1, ther 
fore, were obtained entirely from experiments wit! 
classes. 

The crystals formed were always very small, often 
barely visible under the microscope, so that it was 
usually difficult to determine the refractive index 
even approximately. In particularly 
with the low-silica glasses, hydration and erystalliza- 
tion appeared to take place without altering the 


some Causes, 


shape of the glass fragments 
Table l presents the results of two series of experi- 
ments A) those in which the Morev-tvype bombs 


( AS. 


oe 








sed, and B) 
is Was used. 
order of increasing 
ising temperature, 

Figure 1 is based on the data in table 1, 
mission of experiments in which 
vidence that equilibrium was not reached 
the unit cell size (a function of composition Is 
plotted against temperature of formation. The rela- 
tion between composition and unit cell size is prob- 
ably not quite linear. The composition lines adjacent 
to the vertical seale of figure 1 are based on Yoder’s 
diagram [3], to which reference was made earlier 

Attainment of equilibrium was slow. For exam- 
ple, in experiment 9, after 13 days at 350° C, the 
X-ray pattern indicated the two dis- 
tinet hvdrogarnet compositions, with cell constants 
12.42 and 12.30 A, respectively. This is evidence of 
nonequilibrium conditions. With the higher pres- 
sures attainable in the cold-seal apparatus, equilib- 
rium was reached more quickly, but even here there 
was evidence of variable composition, in that the 
erystals often contained nuclei of lower refractive 
index than that of the surrounding shell This 
indicates a relative deficiency in silica in the nuelei 
of the ervstals, from which it is assumed that these 
may have formed at lower temperatures, during the 
warmup period. In most the amount 
insufficient to affeet the X-ray diffraction pattern, 


those in which the cold-seal 
Within each series the arrange- 
s in silica content and 
with the 
there is 


those 


presence of 


cusses 


Was 


but it was often enough to interfere with a precise 
determination of index of refraction For this 
reason, and also because an extraneous siliceous 


phase usually Was present as a product of the reac- 
little success Was achieved in correlating thre 
index of refraction with the COMposition derived 
from the X-ray pattern 

Considering first the data for the 
moles of silica that is, having the COMM position of 
grossularite it is seen that at the two highest tem- 
peratures a product was obtained having a unit cell 
11.85 A The corresponding figure reported 
for grossularite is 11.84 A 1], or LL.S5 [3] The 
pattern for experiment 41 was prepared with the 
inclusion of an internal standard (tungsten 
the lattice constant given for this product Is probably 


tion, 


with 3 


class 


edge ol 


hence 


L - * 
FiGgure 1 Relation between te mperature of formation and ut 
dae iD hudrogqarnets 7 fhe systlen (‘ald ALO) 
SQ) H.O 
I ‘ er ( ,_ s H «det ‘ ( ) ALO ( Hi.) 





the ( 
fineness of the material it was impossible to 
satisfactory determination of refractive index. The 
loss on ignition was zero, within the limits of analyt. 


correct to +0.01 A. Because of treme 


I iake au 


cal error. This fact, considered with the X-ray 
data, makes it evident that the product was ver 
close to the composition of erossularite ; 

At lower temperatures of formation the cel| COn- 
stant is progressively larger, indicating that the com- 
position is poorer in silica and richer in water 
Here the X-ray patterns are weaker, the lines Often 
broader, and lines of other crystalline phases appeay 
together with those of the hvdrogarnet This js 
understandable, as the silica content of the hydro. 
garnet is now less than that of the glass. The excess 
silica reacts with part of the lime, thus reducing thy 
amount of the latter available for hvdrogarnet 
formation. Some alumina is thereby left unae- 
counted for, but the form in which it is present is not 
known. The identity of the calcium silicate may by 
considered established in at least a few cases. Thy 
product of experiment 51 gave an X-ray pattern of 
moderate intensity, which, after subtraction of th 
hvdrogarnet lines, agrees rather well with that 
given by Heller and Taylor |S], for CSH TF (CSH(B 
according to Bogue’s classification Hy) With this as 
il standard for comparison, the presence ol Much 
smaller amounts of the same phase was detected i 
several other preparations. Ina few cases, however 
the extraneous lines did not fit this pattern 

At lower temperatures — the 
progressively slower. Thus after 2 days at 195° ( 
experiment 30) there was little evidence of altera- 
tion of the glass Glasses poorer in silica were found 
to react more readily, and were therefore used in thy 
experiments at the lower temperatures. The siliea 
content of the hvdrogarnet uppeared to be inde- 
pendent of that of the parent vlass compare, for 
example, experiments 29 and 31 except that the 
hvdrogarnet Was never richer Ith silica than the class 
This limitation accounts for the points in figure 


reaction becomes 


that fall significantly above the curve. Such 
departure should not oecur if the mixtures hav 
reached equilibrium (assuming that the curve is: 


true equilibrium line Probably the explanatior 
will have to await a more complete study of the 
quaternary 

Even with the low-silica glasses, 
lowest temperatures rathes 
The pom! plotted in figure | for 100° © 
hvdrogarnet prepared not from a glass but from a 
boiling mixture of lime water, alumina, and sodiun 
silicate, with the silica in moderate excess 

The use of the low-silica glasses at the higher tem- 


svstem 
the results at th 
were unsatisfactory 


represents a 


peratures is also limited by the breakdown of th 
hvdrogarnet phase to 4CaO-3 ALO,-3HLO (see expert 
ments 2, 3, 6, 10, and 13 

To determine the effect of varving pressure on th 
composition of the hvdrogarnet, experiments 33 ane 
35 were conducted at lower pressures than the other 


convert t ising the lette { {. s,andt 
nd HeO, respect \ followed througt 
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memb rs of this series. The pressures in the Morey 


ymbs (experiments | to 20) were likewise lower. 
in figure 1, the experiments at 1,020 atm are repre- 
ented by filled symbols; those at lower pressures by 
pen & mbols. The fact that a single curve approxi- 
nates both sets of points indicates that pressure is 
frelatively little effect within the range covered. <A 


jmilar conclusion was reached by Yoder on the 

asis Of results at higher temperatures [3] 

3,2. Hydrothermal Reactions of Heterogeneous Lime- 
Alumina-Silica Mixtures 


In the experiments deseribed below, the solid 
omponents were present initially not as a homo- 
veneous glass but as mixtures of various compounds 
In most Cases the results departed more or less from 
hose of the comparable experiments with glasses, 
is described This is aseribed to lack of 
ittainment of equilibrium in the heterogeneous mix- 
ires The results obtained, nevertheless, may be of 
nterest, as the experimental conditions are more 
losely related to that may eXist in actual 
ractical processes 

The data presented in) tables 2? to 5 are represent- 
tive of a considerably larger number of experi- 
ents The data riven were selected hot only to 
llustrate veneral trends but to point out anomalies 


abov °°. 


those 


is well 


Mixtures of Tricalcium Aluminate and Tricalcium Silicate 


The hydration of aluminate in the 
presence of tricalcium interest 
ecause it represents, in greatly simplified form, the 
eactions of the aluminate phase of cement clinket 
luring hydration. In_ the experiments 
ecorded in table 2. the ratio of C.S to C.A ranged 
The ratio is well below that eXisting 
n normal portland cements, but the present study 
sconcerned primarily with the hyvdrogarnet phase, 


tricalelum 


silicate is of some 


SeCTLCS of 


from 0.7 to 3.0 


nwhich the ratio of silica to alumina does not exceed 
1} 

Several facts may be Perhaps the most 
bvious is that in most cases equilibrium was not 
ittained even after 35 experiments 4, 11 
This Is ey ilenced by the coexistence ol more 
vdrogarnet phases. Furthermore, a comparison of 
xperiments Nos. 3 and 4 
toward equilibrium between 1 and 
the experiments helow 170° ¢ and 
is IS2° C (experiments 15, 21) there was evidence 
f only one hvdrogarnet prise and this had the unit 
ell size (12.57 A) of tricalcium aluminate hexahy- 
irate. It must be borne in mind that the presence 
r absence of the hvdrogarnet phase was 


noted 


davs 
two ol 


indicates no 
35 davs 


progress 
In all 


n some as high 


second o 
letermined by examination of the X-ray patterns, 
and is au result the presence ol such il phase in small 
mounts would fail to be Nevertheless, 
t appears that the aluminate and _ silicate phases 
ave hivdrated separately at the lower temperatures 
On the other hand, 
slow as 201° C 
only the siliceous hvdrogarnet Is 
1 dav experunent l7 


observ ed 


above 225° C, and in some cases 
3C'aO-ALO;-6H.O is absent, and 
observed This 
may occur even after only 


The unit cell size of the siliceous hvdrogarnet 
varies little with the quantity of silica in the reaction 
mixture, and apparently not at all with temperature, 
except for a few random On the 
basis of Yoder’s diagram, previousl\ cited, the 12.40 A 
unit cell would indicate about 0.5 mole of SiO, pet 
mole of 3CaQ-AlO;, and 12.32 A would correspond 
to about O.S mole of SiO,. In no case would this 
account for all the silica present; hence part of the 
silica must be present in another phase. In two 
instances (experiments 9 and 22) there is evidence of 
a small amount of hydrogarnet with somewhat 
higher silica. but these appear to be exceptions 
Experiment 8 also fails to conform to the general 
pattern 


mconsistencies 


rarnteE 2 Hlyudrothe mal reactions of mixtures of calc ; 
aluminate and tricalcium silicate 
- 
Phases present pr i 
« . Hvar r 
—_ = — S i] 
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i be CSH (CCSHICA rding to B 1 
Ces hydrate (CeSH( rding to B 
Fron ‘ tenes nyar 


All the reaction products contained free caleium 
hydroxide, us would be expected In most cases the 
X-ray patterns showed a few weak or medium lines 
that could not be attributed either to a hvdrogarnet 
or to Ca(OH), Unaccountably, such lines were 
conspicuously absent from the patterns of experi- 
ments 7 to 12, in which equimolar amounts of the 
aluminate and silicate were initially present, although 
with both higher and lower proportions of silica the 
lines were observed. Interpretation of the patterns 
Is highly speculative, because most of the lines 
were very weak, and undoubtedly many still weaker 
lines were masked by the patterns of the predominant 
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lt is reasonable to assume that they repre- 

or more hvdrated calcium silicates. Most 
extraneous patterns were of two types, here 
lesignated “a’’ and “b.”’ The former bears some 
resemblance to that of the CaO-SiO,-H.O deseribed 


by MeMurdie and Flint {10}, designated CSH (A) 
by Bogue [9], although the fit is imperfect. The 
“b” pattern resembles that of the so-called C.S 


y-hvdrate, which Bogue calls ( SSH 6 Here agath 
the fit is not perfect, but better than in the previous 
The possibility remains, of course, that these 
patterns are indicative of one or more quaternary 
compounds or solid solutions of unknown composition. 
Although the identity of the calcium. silicate 
hvdrates formed in reactions cannot be con- 
sidered well established, it is significant that there 
is no evidence of the presence of either xonotlite or 
hillebrandite, which Peppler 1] has shown be 
the stable phases in the ¢ ‘aO-s1f ).-H () svstem under 
equivalent conditions of temperature and pressure 


Cuse 


these 


to 


b. Mixtures of 3 CaO. Al.O,.6H.O with Silica or Calcium 
Silicate Hydrate 


It is of interest to ascertain how tricalcium alu- 
minate hexahvdrate, prepared in the of 
silica, reacts on subsequent treatment in the presence 
of silica. The data presented in table 3 are typical 
of the results of a considerable number of experi- 
this nature. In experiments | to 10 the 
silica was present as silica gel. As previously noted 


by Flint and Wells [2], the aluminate takes up silica 


absence 


ments ol 


under these conditions, forming silica-bearing hv- 
drogarnets Flint and Wells relied on refractive 
index measurements in estimating the amount of 
silica taken up. The X-ray diffraction method 
emploved in the present work provides 2 somewhat 
more precise indication of the silica content It is 
clear that experiments 3, 5, and 6 failed to reach 

PARLE 3 Hydrothermal reactions | CaO 


/ 





equilibrium, because a new hydrogarnet 


. p USE jp 
seen to exist together with part of the origing 


3C'aO- Al, O,-6H,0. By the use of higher emper. 
atures and pressures, equilibrium apparently was 
obtained. The product of experiment 10 was of 
almost the same composition as that obtained by 
treatment of the ternary glass under similar cond). 
tions. As in the experiments with the 
several of the reaction products formed at the lowe 
temperatures contained hydrated calcium. silieat, 


A lasses 


CSH(B). In a few instances the extraneous lines 
in the X-ray patterns could not be positively 
identified 

In experiments ll to 13 the sihea was initially 


present as a synthetic gel of approximate compositio 
CSH, whereas for Nos. 14 to 16 a synthetic CLSH(B 
containing about 10 percent of xonotlite was used 
These preparations reacted readily with the alumi- 
nate, forming Ca(OH), and a siliceous hvdrogarne; 
No consistent relation is observed between COMDPosi- 
tion of the hydrogarnet and temperature of forma. 
tion; also there is evidence of failure to reach equ- 
librium, 


c. Mixtures of Calcium Hydroxide and Dehydrated Kaolin 


Among the possible practical aspects of hydro. 
formation is the reaction between lime ani 
particularly if conducted a 


garnet 
pozzolanic materials, 
elevated temperatures.  Pozzolans usually 
both silica and alumina in active form. A sample of 
Georgia kaolin, activated by heating at 700° C, was 
used as a representative material of this tvpe. Tabk 
t vives the results of a few experiments with mixtures 
of lime with the activated kaolin The ratio of re. 
actants was such as to give an over-all COTLPOSLUO! 
of 3CaO: ALO,: 2510, 

Like the other mixtures described above, the lime. 
kaolin mixtures failed to reach equilibrium during the 


eontar 
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l period with the possible exception of experl- 
If in the product of each of the last four 
the phase highest in 


react! 


ment 
mem! rs of this series only 


jlica (smallest in unit cell dimensions) be con- 
jdered, the agreement with the equilibrium curve in 
foure is good It should be noted that the unit 
ell size in experiment 7 indicates nearly 3 moles of 
jlica in the hvdrogarnet, which value is higher thar 
the proportion originally present. The excess lime 
ynd alumina are accounted for by the Ca(OH 
nd 4CaO-3ALO,-3H.O present 
TABLI } Hlyudrothe ma reactions oft mirtures oft wie and 
H 
“Hl Ht 
’ ; CSHU 
i 
. . LH 
d. Mixtures of the Hydrated Oxides 
In another short series of experiments, the = re- 
etants were present as hverated oxides: the lime 


is hydroxide, the alumina iis cibbsite, and the silica 
sgel. The ratio of lime to alumina was kept at 3 to 
as in the foregoing, but the silica was varied. The 
Lack of attainment 
reaction period is aga 
The 

3H.LO in experiments | and 3 
It will be shown later that this phase 


esults are presented in table 5 
f equilibrium during the 
bserved im certain Cases occurrence ol 


(‘'agO- 3ALO 
omunernt 


deser\ cs 


vould be expected toappeal only at somewhat higher 


temperatures in the presence of silica in the amounts 
ndicated in table 5.) This may be taken as further 


vidence of lack of equilibrium 


Hlydrothe 
H 
( i 
CAH 
Os 2.4 CAH 
4 j 24 C\AcH 
‘ 1244 Cy,AcH 
é é © i; aad ’ 
! “ is 
k ie i 


3.3. Hydrothermal Decomposition of the Hydro- 


garnets 


As noted in the introduction, trical im 
hexahvdrate decomposes unde! hvdrot! ermal cor 
ditions at about 215° C, vielding calcium hvdroxid: 


and i less basic aluminate ‘| he equation mav be 
written 
3(8CaO-ALO,-6H.O) @4CaO0-3 ALO,-3 HO 

5Ca(OH LOHLO 


Data obtained in the present work (experiments 
and 2 in table 6) indicate that the equilibrium 
temperature is between 220° and 226° C. If a 
hvdrogarnet) of composition 3CaQ-ALO,-mSiO 
§—2m)H,O breaks down in the same manner the 
silica liberated may be expected to combine with 
part of the lime 

Table 6 presents data obtained on hvdrogarnets 
compositions covering a wide range of 
silica content The hvdrogarnets experiments 5 to 
11) were prepared from boiling solutions as described 
in section 2.3. The. silica indicated was 
determined by analysis of the product, but it is 
probable that not all the silica was actually in the 
hvdrogarnet phase, especially in the case of the most 
highly siliceous preparations (experiments 9 to 11 

\ll the preparations, with silica ranging up to 2 
moles per mole of CyA, when treated hydrothermally 
at sufficiently high temperatures vielded decomposi- 
tion products in which (aO-3ALO 3H () could be 


glass 


and 


content 
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identified. Caletum hydroxide was also observed 
except in the most. silica-rich preparation, where 
presumably the lime had all combined with the silica 
One or more other phases, usually ill-defined, were 
also present in the preparations These 
presumably were calcium silicate hydrates, but in 
most cases thes could not be definitely identified 
The minimum temperature of formation = of 
{faQ-3ALO,3H.O is shown as a function of silica 
content in figure 2. The filled circles, representing 
data for the class COM positions, fall close tou smooth 
curve; nevertheless, the diagram should be inter- 
preted with eaution., because complete decomposition 
Was not obtained except ut considerably higher 
temperatures. Furthermore, the open circles repre- 
senting the data for hvdrogarnet starting materials 
are considerably above the curve, whereas under 
equilibrium conditions they should coincide. As a 
possible explanation for the discrepancy, it Is sug- 
gested that the hvdrogarnet phase formed from the 
glass during the warmup period might have been 
poorer in silica than that formed at the equilibrium 
temperature, and therefore subject to decomposition 
at a lower temperature. If this supposition is cor- 
rect, the open circles more closely represent the true 


siliceous 


equilibrium curve 


3.4. Experiments With Cements 


In view of the possible relation between these reac- 
tions and the setting of hydraulic cements, a few 
tests were made on two cements of domestic manu- 
facture, one aluminous and one portland, to see 
whether the hydrated aluminate phases could be 
identified. The anhydrous cement was mixed with 
water and hydrothermally treated for 1 day 





After treatment at 185° C the aluminous ement 
gave an X-ray pattern in which the lines of a hydgp. 
garnet with unit cell about 12.55 A could be dig 
tinguished. Within the limits of measurement, this 
could be interpreted as indicating pure tricaleiyy 
aluminate hexahydrate or a hydrogarnet containing 
a very small amount of silica. The other r 
the pattern fit the published patterns of afwilljy, 
fairly well, but were too weak for positive identifies. 
tion A companion sample treated at 250 
a weaker hvdrogarnet pattern plus a 
lines due to 4CaO-3ALO,.-3H.O The lattes Phas. 
was also observed under the MICroscope Th 
X-ray lines attributed tentatively to afwillite » 
mained the same 


hes “i 


( ; showed 


It is evident that the hydration of the aluminoys 
cement at 185° C produces a hydrogarnet containing 


little or ho silica, despite the presence of sufficient 


silica in the cement to form a hvdrated silicate lt 


may be noted that similar results were obtaine 
from the mixtures of tricalcium aluminate an 
tricalcium silicate at comparable temperatures (se 
for example, experiments 15 and 21 in table 2). Ty 
similarity, however, seems to have disappeared 


150° C 


The sample of portland cement treated hivdro- 


thermally at 185° C gave a product showing a very 
weak hyvdrogarnet pattern with unit cell 12.41 4 
indicating the presence of silica in the hvdrogarnet 
phase. Calcium hydroxide was present in large; 
amount, and the rest of the lines resembled the 
published patterns for 3CaO-SiO,-2H,0. ‘The samp 
held at 250 
tion of formation of 4CaQO-3ALO,.-3H.O. Thus thy 
portland cement differs from the aluminous cement 
in respect to hydrothermal hydration of its aluminat 
constituent It does not necessarily follow that al 
portland cements would behave in similar manner 
The cement contained 7.1 percent of ALO 
which is higher than average. With lower alumine 
content it might prove impossible to identify any 
hydrated aluminate by X-ray means 

The results obtained with both cements should by 
interpreted with caution because of the undetermine 
effect of the iron, alkalies, 
present mm the cement Substitution of FeO for 
ALO 
might 
eontent 


used 


compensate for the effect of higher silica 


4. Summary 


Hydrothermal treatment of lime-alumina-silic 
vlasses of the composition of grossularite produce: 
hvdrogarnets increasing in silica content with th 
temperature of formation and approximating th 
composition of crossularite at SOO? C The curv 
was extended at the lower temperature end throug! 
the use of glasses poorer in silica. 

Mixtures of lime, alumina, and silica as such ort 
various states of combination, on hvdrotherma 
treatment, usually failed to reach the equilibriun 
state attained with the glasses. At times as man} 
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number of 


C gave a similar pattern, with no indica- 


and other components 


would increase the unit cell size [1] and thus 
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distinct hydrogarnet were 


compositions 
to coexist after reaction periods of several 


gs three 


foun 
davs 
TI equilibrium temperature for the decomposition 
of 34 1O-ALO 6H.O undet hydrothermal conditions 
was 29() and 226° C, the prod- 
eing 4CaO-3ALO,-3H,0 and Ca(OH ),. Hydro- 
amounts of silica 
decomposition at 


| sund Lo be between 
cts 
mcreasmg 


garnets containing 
type ol 


mderwent the same 
progressively higher temperatures 

(Co npounds of the hvdrogarnet series were 
d by hydrothermal treatment of aluminous and 
portland cements The product of the aluminous 
ement contained little or no silica, that 
yf thre portland cement contained an appreciable 


amount 


also 


form 


W hereas 
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Infrared Absorption Spectrum of Trimethylborane 


James E. Stewart 


Lhe l ire] iosorptlo 





<pect ‘ guscous ! i BCH 
»ibserved the 2- to 40-micro spect il egio \ nti il issigl ‘ ts i t l 
for infrared bands is well as for the Raman lines reported by othe nvestigators Phe 
perfect wus ¢ tropy Vas eule iluted for i tem rature of 190 G2 kK ind compared Wit t 
experi menta iline Barriers to internal rotation of methvl groups of about 750 « 1.640 
culories pe Me ure culeculated for st ct with or without a threefold axis of s t 


1. Introduction 


study of the infrared 


This 


trimethvlborane was undertaken as part ol a 


spectrum Ob gaseous 
pro- 
gram of the National Bureau of Standards for thy 
nvestigation of chemical and phvsical properties ol 

poron compounds 
The infrared and Raman 
borane were first reported by Goubeau 
» 


trimethvl- 
and Becher 
Their data were used by Siebert [2] in a normal 
treatment. The infrared spectrum. of 
Goubeau and Becher does not extend bevond about 
and it unpurity bands not 
served in the present study \ccordingly, 
desirable to repeat and extend the measurements 
Verv recently Levy [3 the 
nfrared spectrum of  trimethviborane 


horon compounds 


spectra ol 


eoordinate 


H0U em contains ob- 


it seemed 
measured neal 
othe 


very 
and 


2. Experimental Procedure 


TI t suimples of trim thvlborane were prepared by 
\. R and G. oS the 
Bureau's Pure Substances Section by reaction of 
hydrogen = chloride ammonia adduct 
BIC): NH \fter subsequent purification [4], 
purities of two lots were determined calorimetrically 
by G. T. Furukawa of the Thermodynamics Section 
and found to be 99.7+0.1 and 99.9+0.1 percent 
These lots were indistinguishable spectroscopically 

The material was stored as liquid at reduced tem- 
perature in stainless-steel evlinders and was distilled 
into an auNiliary glass bulb immediately before ad- 
mission into the absorption cell. Purification pro- 
cedures were monitored by obtaining spectra after 
Operations on samples from both the 
quid and gas phases In the evlinder This pro- 
cedure was of value in the identification of impurity 
bands mereur\ 
manometer 

Spectra were obtained with a Perkin-Elmer model 
21 double-beam spectrophotometer equipped wit! 


Glasgow. Jr Ross ‘. ol 


Wit! the 


SUCCESSIVE 


Pressures were measured with a 


of sodium bromide for 
the 2- to 15-u and 15- to 40-u regions, respectively 
stray 
bromide prism was limited to 2 
wavelengths less than 3S u by means of reststrahlen 
filters of lithium fluoride and ealeium fluoride \ 
Perkin-Elmer model 12B spectrometer with a lithium 
fluoride prism ised for the near infrared 
region \ potassium bromide prism was not. re 
quired because no bands were observed between 15 
and 25 , Standard Perkin-Elmer gas-cell bodies o! 
nickel-plated brass were used. Potassium bromide 
and bromide were mounted with 
silicone rubber gaskets coated with silicone grease 


prisms chloride and cesium 


radiation encountered when using the cesium 


pereent ol less for 


Wis also 


cesium windows 


with polvchlorotrifluoroethvlene Wax, Ol with a 
beeswaX-rosin mixture 
3. Results 
Trimethylborane Is verv reactive with oxvgen 


and ignites spontaneously in air. Very slow decom 
position of the sample in the cell was observed, 
probably due to slight air leakage at the window seals 
No spectral evidence was found for the formation 
of methane, ethane, ethylene, or dioxide 
during this The appearance of a band 
near 1370 em! is very likely associated with the 
formation of B—O—C bonds [5], suggesting that 
boron methonides are produced by slow oxidation 
Additional bands appeared at 1,479, 1,220, 1,062 
1.055-1,048 (PQR-like structure), 950, 787, and 590 
cm When purified material stored for 
several weeks at reduced temperature and reexam- 


Cul bon 
reaction 


Wis 


ined, none of the above bands was seen exce pt 787 
em An additional band was found at 730 em 
Goubeau and Becher reported medium to strong 


620. 635. 648. 667. 720. and SOO em 

which were not observed in the present study even 
at higher Furthermore, 
were not resolved in their spectrum because suffi- 


crently low pressure Was not used. 


bands at 


pressures several bands 


The sodium chloride and cesium bromide spectra 
are shown in figure 1, and the complete spectrum 
is tabulated in table 1, along with the Raman spee- 
trum of Goubeau and Becher. 
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C,B bendit 
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B 
oCH ( 
oCH i 
O6SN P brane} 
CH voli 
O75 320 Quay 
lmaip 
CH 
S55 so 1173 
S55 rot 110] 
Y8S3 + 320 = 1503 
UOHS ul 304 
KC 
CH t 
(‘Hi i 
2X 85d 1710 
1149+ 675 IS24 
13090 +675 1OS4 
1300 S55 -1o4 
Al ee | O75 2 4H 
1459 + 983 — 2442 
1300-4 L149 = 245s 
1162 1300 — 2471 
1149 1440 —- 2580 
vCH 
vCH 
vCH i 
vCH ‘ 
PQ75 520 = 8295 
2875+ 2 320=— 3513 
AL ee | 2 320 Sob! 
SOD 2875 3730 
968 + 2875 — 3845 
1149+ 2921 1070 
2986 - 1149 $135 
2975 1162 1137 
2921 1440 ie 
2921 + 2875 — 5796 
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4 Structure, Vibrational Assignment, and 
Discussion 


Electron diffraction studies [6] have shown. tri- 
ethvlborane to have a planar C,B_ skeleton of 
ymmetryv J). The configuration of the hvdrogen 
loms is not known Barring completely free 
thvl group rotation, the molecule can be in one of 
¢ following symmetry classes: Cy, with one hydro- 
n from each methyl on one side of the C.B plane 
d the remaining hydrogens on the other side; 
with one hydrogen from each methyl in the C,B 
with the methyls rotated by an equal 
jount about their axes; C, with a plane of symmetry 
ntaiming the ¢ axis of the C.B skeleton: or ¢ 
th no symmetry at all. The study of molecular 


ane; f 


dels shows that the (, strueture can probably 
ruled out, whereas the (Cy, structure Is most 
tractive from steric considerations. The number 
i vibrations observed in the Raman effect but not 


Is @\ idence against 
transitions 


vice versa, 


which all 


the infrared, and 
(., and (. symmetry, in 


are both Raman and infrared active. On the other 
hand, the abundance of vibrations observed in the 
850 to 1,200 ecm region makes the nondegenerate 
(’; or Cy symmetry more reasonable than the de- 
generate structures. In view of these considerations, 
both the (, and the Cy, structures are treated, with 
the reservation that in the case of C, symmetry the 
C.B skeletal modes may behave as though they were 
in a molecule of Jy, symmetry because of weakness 
of methyl-skeleton interactions. 

There are no. assignments 
boron-alkyl compounds to be used as a guide in the 
analysis of the trimethylborane spectrum, Sheppard 
and Simpson |7] have considered the effect of sub- 
stituting fluorme for methyl groups ia several 
paraffins and found good correlation between skeletal 
vibrations with usually little frequeney shift for the 
unsvmmetric vibrations and a moderate shift to 
higher frequencies for the svmmetric vibrations of 
the fluorine analogs. This suggests that the vibra- 
tional assignment for BF, can be used as an aid in 
the skeletal trimethylborane 


available for other 


modes of 


ussigning 
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to demonstrate the extent of the reliability 
procedure, data for five pairs of compounds 
mbling BF, and B(CH and not discussed by 
and Simpson, tabulated in table 2. 
Phe ske letal modes for neetone are assigned from the 
frequencies given by Price. All other assignments 
are those given by the original author. With the 
exception of the asvmmetric stretching modes of 
NF, all vibrations are shifted to higher frequencies 
upon substitution of fluorine, some by quite large 
nmounts Of the two possible choices for the ausVvim- 
metric stretching vibration for trimethylborane at 
1.309 and 1,149 em, 1,309 em seems to be more 
probable by analogy with boron trifluoride. This ts 
in disagreement with Siebert’s choice of 1,149 em 
The infrared spectra of a series of seven trialkvl 
boranes have been measured by Pickard, Dulaney, 
and Polly [S]. In addition, the infrared * and Raman 
9} spectra of triethvlborane have been measured at 


The ypard ure 


the Bureau Krom these curves it Is possible to 
select a satisfactory sequence of frequencies con- 
verging on 1,309) em for trimethyvlborane for 
assigninent to asvimetrie skeletal stretching vibra- 
tions 

Intensity considerations supply further evidence 


concerning the assignment of the asvmmetric stretch- 
ing vibration. The band at 1.309 em~! must include 
some contribution from the svmmetric CH, defor- 
mation modes Sut it much Intense 
relative to the asvmmetric CH, deformation at 1,459 
to be entirely due to CH vibrations Hence 
vibration Further- 
in the Raman 


sSeenis too 
Cin 
another must be 
more the great Intensits of 1,145 em 
effect argues against its assignment to a nontotally 
SVinmetriec mode The corresponding vibration im 
the Raman spectrum of BF, is observed even 
though it forbidden transition. On the 
other hand, no band is available in the neighborhood 


supel posed 


hot 
is hot a 


of 1,309 em~! for assignment to the B" isotopic mole- 
cule, whereas, 1,149 em has such a neighbor of 
proper intensity at 1,189 em 

Siebert’s assignment of 486 em~ to the BC, out- 


of-plane bending vibration is satisfactory by analogy 
with the corresponding frequency in BF, but this 
vibration is Raman inactive and infrared active under 
the DP BF. (and also under (Cy 
whereas just the reverse activity is observed for the 
$86 em”! band. Therefore, em”! is assigned to 
this vibration The assignment of 320 em to the 
in-plane deformation mode is natural and agrees with 
Siebert 

A choice must be 
O75 and S64 em lor assignment to the svmmetric 
B—C stretching mode. 864 em~' is perhaps in better 
agreement with the BF, analogy, but its use would 
leave 675 em~! unexplained, whereas 864 ¢m 
reasonably methy| mode 
Is reported to be sharp ana 
strong in the Raman effect ana Is unobserved 


selection rules of 
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made between the vibrations at 


can be 
assigned to a rocking 
Furthermore, 675 em 
very 
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1€8 OF fluoro / 


TABLI » Nke ela ribration i equen 


methul compounds 


. \svm- Svm- Out-of- | slang 
Compound 
a ee metric metric plane for- 
stretch stretch bending matio 
Cri hii Ahi 
NF 1, 031 QO7 642 197 
N(CH 1, O36 S27 $235 65 
PF RQ? S60 1S7 14 
P(CH TOS 653 305 262 
\sF 707 644 341 274 
AS(CH DS5 56S 239 Lee 
O—CF 1. 249 O65 626 584 
O=—-C(CH “ 1, 202 SOG 47 
F.C=CH 1, 302 926 611 550 
CH.);<C=CH OS6 SOO 30] S78 
BF 1, 446 SSS Ho] iSO 
BiCH 1, 3090 675 336 20) 
Kk. L. Pace and Louis Pieree, J. Chem. Phys. 23, 1248 
1055 
Kk. W. F. Rohlrauseh, Monatsh. 68, 349 (1936 
MI. K. Wilson and S. R. Polo, J. Chem. Phys. 20, 17 
1952 
tk. J. Rosenbaum, D. J. Rubir ind ©. R. Sa rg 
Chem. Phys. 8, 366 (1940 
J. B. Howard and E. B. Wilson, Jr J. Chem. PI J 
O30 (1934 
\. H. Nielsen, T. G. Burke, P. J. H. Wo 1 | 
Jones, J Chem Phvs 20. 56 1952 
D. Price, J. Chem. Phys. 9, 725 (1941 
D. C. Smith, J. R. Nielsen, and H. H. Claassen, J. Cl 
Phys. 18, 326 (1950 
J. E. Kilpatrick and K.S. Pitzer, J. Research NBS 38, ] 
1947) RPI76S 
1) \l Crug ind E I Barker J. Chem Phy 47. 453 
1939); D. M. Yost. D. DeVault, T. I \ndersor 1 | 
Lassetre, J. Chem. Phys. 6, 424 (1938S 
in the infrared, whereas S64 em~! is broad in ti 


Raman effect and of medium intensity in the infrar 
spectrum. The assignment of 675 em 
metric stretch is therefore preferred, in agreemer 
with Siebert 


to the SVN 


The skeletal force constants were calculated I 
treating the methyl groups as point masses [10 
The B & stretching force constant lac. which d 

! 


pends on the asvmmetric stretching frequency onl 
is found to be 4.03 « 10 dynes em as compared wit 
Siebert’'s 3.04 10) dvnes cm The Sell 
consistency of the point-mass methyl treatment ca 
be demonstrated by comparing the caleulated val 


value of 


of (y with the experimentally determin 
value. Using Siebert’s assignment, these number 
are 1.57 10° and 1.42 10 em respectivell 
The modified ussighniment proposed here o1Ves 
1.6250 10° and 1.82 10°, respectively 


The methyl stretching and deformation modes a 
assigned in the usual way. Assuming a structure | 
( svinmetry, the bands at 855. 906. and 1.149 em 
are assigned to a’ methyl rocking modes, and 6S 
983, and 1,162 em are 
modes These bands are assigned in such a Way Uh 


assigned to oa’’ rocking 
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n 
transitions are most intense 
Assuming a structure of 


possible the a’ transitions are most intense 


man effeet and a’’ 


wher 
ne hk 
» the infrared spectrum 

svinmetry, the methyl rocking modes are as- 


a’’-968, e’-1,149, e’’-864 


violate the selee- 


follows: a’-906, 
m though 968 and S64 em 


is 


on rules. Raman lines at 394, 486. 755. and 1.060 
Mm emain unexplained. The 1,060 em ~! line 
aight correspond to the impurity band observed 
ear this frequency in the infrared spectrum. The 
nes at 486, 755, and 1,060 em”! might arise from a 
vethv! torsion vibration of frequency 160) em 

9) O00 SO. GOGH 160— 746. 906 160 1.066. 
fhe shoulder at 345 em~' in the infrared spectrum 
night result from a splitting of the skeletal deforma- 
on vibrations, or it might represent a B" isotopic 
nolecule vibration. The weak shoulder at about 
60 em! in the infrared spectrum might be the P? 
rane! corresponding to 968 em—', or it might be an 


ipper stage band 96S In the assignment 


or the ¢ structure, the combination transitions 
sumed for the weak bands at 1,300, 2.444, 2.472. 
nd 3.729 cm are in violation of selection rules, 
nd the strong peaks at 9S3 and 1,162 em~! are as- 


to combination transitions 


5. Entropy and Methy] Rotation Barrier 


Using observed vibration frequencies, moments of 
ertia f= In= 1123.92X< 10 and | 209.02 
d-” e-em? and a reduced methyl moment of mertia 
f 5.322 * 10 g-cm’, the contribution to the entropy 
it 199.92° K from translation, over-all rotation, m- 
ernal free rotation, and vibration is 69.3] eal deg 
for the C.», structure and 71.48 eal dee 
nole' for the C, structure. The difference is largels 
lue to the difference in svmmetry number of 81 for 
, and 27 for ¢ Recently Furukawa and Park [11] 
letermined the entropy of trimethylborane calori- 
netrically and found S 68.29 cal deg”! mole 
The difference of 0.34 cal deg~' mole”! per rotor for 

and 1.06 cal deg”! mole for (. symmetry re- 
wires for agreement a sinusoidal barrier to free ro- 
ation of 750 eal/mole or 1,640 cal/mole for the Cy, 
ind (. structures, respec tively 


nole 


6. Conclusion 


In several respects the assignment for a= trim- 
thv borane structure of C. symmetry ts more satts- 
factory than the assignment for a C', structure 
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Furthermore, the higher methyl rotation barrier 
associated with the C, structure would seem more ac- 
ceptable upon comparison with similar molecules, and 
indeed French and Rasmussen [12] predicted a bar- 
rier of 1,720 cal/mole on empirical grounds. On the 
other hand, the nature of the barrier rotation 
about a single bond and the role of methyl inter- 
actions in determining the barrier are not understood. 
Hence it is by no means clear that the use of standard 
methods for calculating the barrier is justified in the 
case of trimethyvlborane, or that the comparison of 
barriers between molecules is meaningful. There- 
fore, made regarding the correct 
structure of trimethyvlborane 


to 


ho conclusion Is 


The author is grateful to A. R. Glasgow, Jr., and 
G. S. Ross, who prepared and purified the sam- 
ples used and performed the cell filling operations; 
to G. 'T. Furukawa and R. P. Park for permitting the 
use of their thermodynamic data in advance of pub- 
lication; to D. Kk. Mann for his helpful discussions 
and advice; W. H. Evans for checking the 
thermodynamic calculations. 
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Samples of cotton. cotton 


a high vacuum The 


of CO, COs, water, and levoglucosan (tar 
lm pregnation 


an increase in the vields of CO, COs, 


same materials were 


Ising a Very sensitive 


of weight 


rials are much greater than those 


or sodium chloride 


l. Introduction 


An investigation of the thermal degradation or 
pyrolysis of cellulose in a vacuum was undertaken 
n connection with the problem of rendering cellu- 
materials used in fabrics flame-resistant. 
Flaming of these materials takes place in two prin- 
ipal consecutive steps | the cellulose first 
grades under the influence of heat to form volatile 
and nonvolatile products; (2) the flammable consti- 
tuents of the products ignite and the heat evolved 
promotes further decomposition of the cellulose. In 
rendering cellulosic materials resistant to flaming, 
the conditions of degradation are so altered as to 
nerease the amount of the nonflammable products 
at the expense of the flammable products. 

This paper deals with the first of the two steps 
numerated above. The advantages of conducting 
pyrolysis in a vacuum, instead of in air or in an 
inert gas, are that the degradation products are pre- 
vented from reacting with oxvgen (in the 
air) and that the volatile products are more quickly 
and completely removed from the hot zone, where 
they might otherwise undergo secondary reactions. 
Ina vacuum the volatile products are more easily 
collected and fractionated for analysis. 

P\ roly SIs of cellulose at reduced pressures has been 
reported by several investigators [1], 2,3, 4 Their 
iidieate the formation of tar, char, H.O, 
CO,, and CO in various vields, depending on condi- 
tions of pyrolysis. In this investigation a vacuum 
£10~ to 107° mm Hg was maintained during pyrol- 
ysis, and a study was made of the nature and pro- 
portion of the various products and also of the rates 
of degradation of pure cellulosic materials with or 
Without additives. A amount of re- 
search has been conducted in recent vears al the Bu- 


oOsic 


de- 


case of 


results 


considerable 





Pyrolysis of Cellulose in a Vacuum’ 


hvdrocellulose, 
mpregnated with sodium carbonate or sodium chloride. 
l volatile products were fractionated and the 
the mass spectrometer and by infrared absorption. The 
The residue consisted mainly of carbon 
of the cellulosic materials with salts caused a decrease in the vield of tar and 
H.O, and char 
investigated in the range 
tungsten spring balance enclosed in a vacuum 
versus percentage of loss of weight, in the case of pure 


through 

ization end point I; the case 

chloride, the initial rates of loss of weight are 
tion end point 


of the same 


> 


S. L. Madorsky, V. E. Hart, and S. Straus 


and viscose ravon, both by themselves and 
were pvrolvzed at 250° to 397° C 
fractions analvzed in 
fractions consisted mainly 


char 


volatile 


Rates of thermal degradation of the 
245° to 305° C bv a loss-of-weight method. 
Plots of rates of loss 


cellulosic materials, pass 


maxima at about 13 to 23 pereent loss of weight, then drop gradually to the earbon- 
of samples impregnated with sodium carbonate or sodium 
very 
The activation energies of thermal degradation of the pure cellulosic mate- 


high, but drop rapidly to the carboniza- 


materials impregnated with sodium carbonate 


reau on the thermal degradation of vinyl, butadiene, 
and aromatic polymers; the methods and apparatus 
developed and used in this research have now been 
applied to the study of thermal degradation of cellu- 
The method can be divided into two parts: 

|) Pyrolysis of a small sample in an evacuated 
apparatus, and subsequent investigation of the de- 
composition products, 

2) Determination of rates of thermal degrada- 
tion at several temperatures, and calculation of the 
activation energies of the degradation reactions in- 
volved 


lose. 


2. Materials Used 


Three forms of cellulose were used: cotton cellu- 
lose, cotton hvdrocellulose, and viscose ravon The 
cotton, which was an Empire cotton obtained from 
the Southern Utilization Branch of the U.S. Depart- 
ment of Agriculture, was dewaxed by extraction with 
boiling alcohol. It was freed from pectic substances 
bv immersion in boiling 1-percent sodium hydroxide 
in the absence of air for 2 hr, washed in distilled 
water, neutralized with dilute acetic acid, again 
washed in water, and dried in air at room temperature. 
The hvdrocellulose was obtained by immersion of a 
portion of the purified cotton in 6 NV HCI for ap- 
proximately 15 hr. After washing with distilled 
water to remove the acid, the hydrocellulose was 
dried in air at room temperature. The viscose rayon 
was a commercial varn, 900 denier, 100 filament. It 
was washed in warm dilute aqueous ammonia, rinsed 
in warm water, and dried in air at room temperature." 

In order to determine the amount of moisture and 
other volatiles adsorbed on the cellulosic materials, 
samples were “dried” to a constant weight in a 
vacuum at 100°C. It was found that this could be 
accomplished in about 60 min. In percent- 
ages of dry samples, were: 3.7, 4.1, and 8.6, for 


| 4OSSCS, 


* The authors are indebted to Florence H. Forziati of the Textiles Section of the 
tor Lose nd rayon 
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and Viscose ravon, 
In the actual pvrolysis expcriments, the 
dried im the manner im the 
lum apparatus and the volatiles were pumped 


The original weights of the 


vdrocellulose., respec- 


above 


~ were 


lo pvrolvsis 


prio 
ipl were corrected for the above determined 
} OSses 
In some of the experiments, samples of cotton 
mpregnated with NaCl or NacCQOs and viscose ravon 
lmpregnated with NVaCO were used The cellu- 
lose Was impregnated with the salt in the following 
manner: A weighted amount of cellulose was placed 
in the contamer in which it was to be pvrolvzed. A 
Weighed amount of salt was added to the same con- 
{aime \ small amount of water was then added and 


the water evaporated in an oven at about SO° C. 
Drving of the sample and additive was then com- 
pleted in a vacuum in the same manner as for the 
All calculations of losses, 
the 


nonimpregnated samples. 


vields residues, etc were made on basis of 


drv celluloses 


3. Pyrolysis 
3.1. Apparatus and Experimental Procedure 


In order to obtain a detailed analysis of the vola- 
tile products of pyrolysis, these were separated into 
fractions. The following volatile and nonvolatile 
fractions were collected: 

|. Residue. 


2. Fraction Vp,,, volatile at the temperature of 
pyrolysis, but not at 25° C 
3. Fraction V.-;, volatile at 25° C, but not at 
Be” ©. 
t. Fraction V volatile at 80° C, but not at 
190° ©, 
5. Fraction Vey, volatile at 190° C 
In previous work [5] on pyrolysis of polymers, fractions Vpyr, Vas, \ 
V_s00 were di , II, LIB, FILA, and LV, respective 
( 
At 
> Sir 1 >) 
-— . He 
| \ 
Kj 
3 ff) 
} | in 
ae - oo va 
P i 0" 0 
4 | " y HW 
3 l 
A /\ 
L 
MM G 
F 
Pieure |, Ipparatus for t ly of the 
A, Small quartz tube in which the samp placed: B, lar 


(a, Ca, tub r collect t tion \ l., trap 


The pyrolysis apparatus is the same as | ised 
in the thermal degradation of polvtetraflucroethy. 


lene [6], and is shown diagrammatically in figure’) 

In figure 1, A is a quartz tube, closed at one end 
and fitted into a larger quartz tube, B, also closed gy 
one end, Tube B is attached to the vacuun sVstem 
by means of a ground joint, C. A vacuum of abou 
10-* mm of mereury could be obtained by means of 
a rotary oil pump, not shown in the figure, » liquid 
nitrogen-cooled trap, P, and a mercury-diffusion 


pump E. A evylindrical electric heater, not showy 
in the fivure, could be moved in a fixed horizontal 
position so that tube A fitted approximately in th, 
center of the heater muffle. The temperature was 
measured by means of a platinum versus platinum-10 
percent rhodium thermocouple fixed permanently jy 
the muffle of the heater, so that when the heater was 
moved in position for pyrolysis, the thermocoupk 
came in contact with tube B approximately at 
point where the sample was located in tube \ near 
its closed end The heater was so constructed that 
an essentially uniform temperature could be main- 
tained over the entire specimen, The temperature 
could be maintamed within 0.5 deg C at pvrolysis 
temperatures . 

After the cellulose specimen was placed in tub 
A and the tube inserted into tube B, the apparatus 
was evacuated overnight. When the heater, which 
had been kept at the desired temperature, but away 
from the apparatus overnight, was moved into 
position, tube A reached the desired temperature 
within a few minutes. 

During pyrolysis the apparatus was cut off from 
the evacuation pump by means of stopcock D, while 
the mereury diffusion pump E remained in opera- 
tion. The combined volatile fractions V... V_.. 
and V were confined between D and EE. Th 
mercury pump was rated to hold 25-mm back pres- 
sure, with a high vacuum on the inlet side: it was 
actually operated ata back pressure of not ove! 
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The pressure was measured in the calibrated 


§ mn 
j space between this mercury pump and stopeock D 
by means of a multiplving manometer, F, containing 
a low vapor-pressure silicone oil in the left arm of 
the manometer, above a mercury column 
A sample of the gaseous fraction \ noncon- 


densable at the temperature of liquid nitrogen, Was 
colle d for analysis at the end of an experiment in 
the following mannet Liquid nitrogen was placed 
trap L, The condensable material, consist- 


aroul 
| ing of combined fractions Veg and V.;, was drawn 
mto L., and a sample of the remaining gaseous frac- 
tion \ was then taken for analysis in 60-ml tube 
G ol (;' by sealing off the tube at IX or K’ The 
welgiil of fraction \ Was calculated from. its 
volume, pressure, and composition 

After evacuating the excess fraction V the 


material that condensed in trap L was separated into 
fractions, V and \. with the aid of dry ice 


two 
| These fractions were collected by means of liquid 
nitrogen in 3-mm tubes M and AI’, which were at- 


tached to the vacuum svstem by means of ground 








joints O and O’. The tubes were sealed at N and 
N’, weighed, and the contents analyzed in the mass 
spectrometer 

Fraction V, volatile the temperature 
pvrolysis but not at room temperature 
during pyrolysis as a vellowish tarry material in the 
quartz tube B just outside the heater. It was dis- 
solved in methanol, transferred into a weighing 
bottle, the solvent evaporated, and the tarrv material 
was then weighed 

The nonvolatile residue remaining in tube A was 
weighed in a weighing bottle. The total recovery 
weight of the products volatilized in each pyrolysis 
experiment was within a few tenths of a milligram of 
the loss in weight of the original drv sample. All 
weighings were made on a semimicro balance 


of 


deposited 


| 


3.2. Results of Pyrolysis 


pyrolysis experiments are shown in 
2 2¢, inclusive, were car- 


ail of 


Results of 
table 1. Experiments 
ried out in consecutive steps on the same specimen 


to 
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y Similarly, experiments 5a to 5f, inclusive, 
arried out on the same specimen of hydro- 
se, and experiments ba to 6f on the same speci- 


These three series of experl- 


mn ol viscose ravon 

ments were carried out at 280° CC. A sample of 
cotton impregnated with Na,CO, (7% was pyro- 
lvzed in two consecutive steps, at 251° and 291° C 
\ sample of cotton impregnated with NaCl (8%) was 
pvrolyvzed in one step at 280° C Fractions V,,, on 
mass-spectrometer analysis, were found to consist 
mainly of water mixed with a small amount of 
nee tale hvde aus 1s shown later Fractions \ and 
\ also on mass-spectrometer analysis, were found 


to consist of CO, and CO, respectively. According 


to table cotton and hvdrocellulose behaved Simii- 


larly on pyrolysis, whereas viscose ravon vielded 
more H.O, CO., and CO, and less tar, than the first 
two cellulosic materials. The addition of NaCl or 
Na-CO, had the effect of mereasing the vields of 
H.O, CO,, and CO, and decreasing the vield of tar. 


In table 2 cumulative vields of fractions are shown 
for 2S0O & except mn the case of purified cotton 
impregnated with NasCOs;, wher the initial rate of 

7% volatili- 


degradation was so fast th; the first 37 / 








- ‘ | | 
zation was carried out 2°251° C and an additional 
14°, at 291° C). Die data shown in table 2 were 
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olan Naot ¢ 
i Al } 2 | ! é a IS oblaine he 
cell 
\ \ \ 
I ( ( ) 
f 
H aN ‘ 2 
| 
Tat Efi f j 


346 





calculated on the basis of the data in table | 


Inter. 


polating where necessary. 
The temperature effect on vields of the irlous 
fractions is shown in table 3, where a comparison jg 


made of results obtained at different temperatures 
but for about the same extent of degradatio: 
According to this table the temperature — effee 
through a ranve of about 250° to 32 C is ne 


pronounced 


4. Analysis of Fractions 


4.1. Mass Spectrometer Analysis 


Mass-spectrometer analysis of fractions \ an 
V _s0 them to CO and CO,, py 
spectively. Similar analysis for fraction — \ 
shown in table 4, indicates that water is the mai 


showed consist of 


TABLE 4 Vass-s pectrometer analysis of action \ 
{ ‘ t 
\ oad te 
{ ¢) 
{ 
Nal ‘ 
VJ VJ VJ VJ 
Wa “4 y 
Ac . 
An o ! “ 
e \ \ 
constituent Che only other constituent is ace. 


taldehvde, which wopears mb amounts varving from 
3 to 16 percent The reason for this variation is not 
apparent from the present limited data 

Analysis of fraction V,,, (tar) was performed b 
Bradt and Mohler in a 60° Nier mass 


method thar thes have developed ri 


sper trometer 


using 2 


The sample was heated in a vacuum-tube furnac 
that opened directly mito the lroOniIizRtIO! chamber 
and the temperature Was increased step by step 
until there was ne further evaporation Mass 


spectra were recorded ut each step Ih) addition 10 
fraction V..., these investigators also vaporized o 
pyrolyvzed directly into the mass spectrometer filt 
paper, purified cotton, viscose ravon. cellobrose. a 
levoglucosan. There 


of pure compounds available for mass-spectromet 


ure ho baekoround patterns 


analysis of these substances, and the conclusions 
drawn on the basis of the mass spectra obtamer 
must therefore be considered as tentative 

Filter papel cotton and viscose ravol Terk 
similar mass spectra When heated above 100° ( 
there was some evolution of HO. CQO., and CO 
At 240° and higher, a rather complicated mass 
spectrum extending to mass 144 appeared. Ce 
lobiose appeared to degrade at 250° CC. and the 
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trlous 


SO] | 
Lures 
LatioT 
eff "1 


is Hol 





spectrum showed smal! peaks, indicating Hiwsses 
eavicer than 144 Its spectrum was much like 
hose f the cellulosic materials 

Le lucosan evaporated in a narrow temperature 


eat about 130° C Its mass spectrum showed 


eaks corresponding to HO, CO,, and CO. These 
eaks were rather small as « ompared with those that 
ere tained for the celluloses 

The mass spectra obtamed for fraction V, tar 
rom tton and = viscose ravon were similar At 

evaporation temperature of about 130° C, these 
spectra were almost identieal with that from 
voglucosan. However, there was in the tar a 


mpo mud that was less volatile than levoglucosan 
nd thiat degraded at higher temperatures to give 
spectrum more like the product obtained from 
irect pyrolysis of cellulose into the mass spectrom- 
tel lt may be significant that all the spectra, 
hethet from. cellulose, levoglucosan, or tar, Con- 
ined peaks for masses 126 and = 144 These 


jasses correspond to that of levoglucosan or of the 
lulose structural unit, molecular weight 162. 
th 2 and | molecules of water, respectively, split 


Infrared Analysis of Fraction V,,, 


Samples of pure levoglucosan and of fraction V, 
Thin 


ms were de posited on rock salt Hy evaporating the 


ar) from cotton were dissolved in methanol. 


lvent vnd infrared absorption spectra of the 
sulting films were obtained As shown tn figure 2,° 
only significant difference between the two 


ectra ois thre presence ol an absorption band | 


4.3. Average Molecular Weight of Fraction V, 


The molecular weight of fraction \ 
tar) from cotton was determined by a microfreezing 
point lowering method from a solution in phenol. 
Four determinations were made on this 
and four on pure levoglucosan as a 
molecular  freezing-point 


wverage 


fraction 
The 
phenol IS 


control 
lowering for 


viven in the literature as 7.3° C Sufficient sample 
was used in each determination to give a freezing 
point depression of approximately 3° C. The 


molecular 
166 2? as 


levoglucosan was found to be 
compared with 162 theoretical. The 
average molecular weight of fraction VV... was found 


to be 152 | 


weight ol 


5. Rates of Degration and Activation 
Energies 


5.1. 


The rates of 
materials were 


Rate Measurements in a Spring Balance 


thermal degradation of cellulosic 
determined by means of a very sensi- 
tive tungsten spring balance enclosed ina Pyrex 
housing, which could be evacuated lo about 10 Inm 
Hg. This apparatus and the method have been de- 
scribed previously [6] in connection with a study of 
the thermal degradation of a number of polymers 
The method, briefly, heating a sample ol 
about 4 to 6 mg in a crucible suspended from the 


COnSIStS In 


tungsten spring, and recording at regular intervals 
the loss in weight of the sample due to pyrolysis and 
evaporation, as indicated by a change in position of 


wire extension of the 
In preparing for a rate 


i crossline Ol a spring 


experiment the crucible 





,78 microns in the speetrum of fraction \ tar , 
| = containing the sample was suspended from the 
is suggests that fraction \ contained both . 
| tungsten spring into a tube projecting down from 
roglucosan and a compound that contains ¢ () ah 
the vacuum apparatus. The apparatus was pumped 
— to a vacuum of approximately 107° mm Hg and left 
| . x 
standing In the meantime an electric furnace, kept 
in a& position away from the projected tube, was 
£ R I heated overnight to the temperature required for the 
B . 
given experiment. In the morning the furnace was 
. 
| 
| SAN 2 dl ~ 
2° 
x | 
= | +FRACTION, Vpyr 
8 
= 20 | 
= . 
WA TH, MICRON 
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raised into position around the sample It ordin 


took nbout 1D) min tor thre 
required temperature 


crucible to re: 


Duration of the experime 


was counted from the time this Lemperalure 


re ached 


The te mperature Was « ontrolled by 


ned 


r ( 


of an electronic thermostat to within ().? «deg 
The results of rate experiments on 
shown in figures 3 and 4 The rates, hk 


age of the original dried cellulose pel 


initially low and pass through 
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wy s heate ) va a ; 0 iHlov ” equated } 
Naol OD) ; nope entagde ( 1108 ( lie ( mie ‘ 
( ono pie entagde Hwa atior A ate 


hich they follow straight lines The pots at 


shich extrapolations of these lines intersect with the 


bscissas are defined as carbonization end points 
xypressed in percentage of volatilization. The char 
eld, in percent of sample, for any particular tem- 
erature can be obtained by subtracting the car 
wnization end point from 100.) By extrapolating 
e straight parts of the Ky rate curves to. the 
dinate, the apparent initial rates were obtained 
e importance of these apparent initial rates lies in 
¢ fact that the samples have the same composition 
ly at zero degradation, and these rates can thus 
used for the purpose of comparison or for cal 
dating activation energies 
The curves obtained by plotting percentage of 
mpl volatilized and logaritl m of percentage resi- 
we against time for hvdrocellulose and viscose 
von, Were similar to the curves obtained for cot- 
hvure 2, nd are not shown here The Ky rate 
fvolatilization curves for hvdrocellulose and viscose 
von are shown tn figures 5 and 6, respectively The 
parent initial AY rates and the carbonization end 
mts for hvdrocellulose and viscose ravon were ob- 
uned by extrapolation us in the case of cotton 
Results ol rate experiments on cotton lniprey- 
ted with NaoCO, (9 ©) are shown in figures 7 and 
The thermal stability ol the cotton was creatly 
duced as a result of the admixture with sodium car- 
Onate This reduetion in stabilits Was particular, 


treat for the bittal stages ol the degradation In 


e lo-min interval during which the cotton was 
used from. room temperature to the desired tem- 


perature, about IS to 30 pereent of the cotton vol- 
atilized, depending upon the final temperature 

The rate curves for Na,CQO,-impregnated cotton 
are shown In figure Ss In this case the rates were ex 
pressed in percentages of residues existing at thy 
Instant (excluding NasCO.) instead of percentages of 
the original samples and are designated as AY rates 
This was done for the reason that the Ky plots are 
straight lines bevond 30 to 40 percent volatilization 
while the A, plots are curved through their whole 
extent Extrapolated initial rates are the same 
whether they are obtained from A, or A, rate curves 

The plots showing percentage volatilization versus 
time and logarithm of percentage residue versus 
time for viscose ravon impregnated with Na.CO 
(96%) and cotton impregnated with NaCl (99%) are 
not shown here. These plots were similar to the 
corresponding plots for cotton impregnated with 
NaCO, (9%), shown in figure 7. The Ay rate 
curves for viscose ravon impregnated with Na.CO 
and cotton impregnated with NaCl are shown in 
figures 9 and 10, respectively. The straight parts 
of the A, plots in figures 9 and 10 were extrapolated 
to obtain the apparent initial rates and the carboni- 
zation end points. —Inall cases, the plots of logarithm 
of residue versus percentage volatilization (shown tm 
figs. 3 and 7 as dotted lines), were not straight lines, 
indicating that the over-all degradation reaction is 
not of first order 
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Taste 5 Rates of thermal deqradation of ce ilo imple 


Pertinent data on all the rate experiments by tl 
loss-of-weight method in the tungsten-spring appar 


tus, are shown in table » Phe activation chergies 
are shown in the last column of this table Thes 


activation energies were obtained graphically 
shown in figure 11, by plotting the log, of the upp 
ent initial rate at any given temperature vers 
1/7 10°. where 7 1s the absolute temperatul 
The over-all activation energy, /, in kilocalories p 
mole. is the slope of the straight line, obtained fror 
the formula 


d log I 


f / | / 


\ comparison of initial rates for all the materials 


investigated Is given in table 6 for 275° and BSS? | 
the two temperatures at which the most complet 
data were obtained Table 7 CIVES il condenst 
summary of data showing the effeet of additions 
NaCO; and NaCl on the vield of flammable at 
nonflammable fractions in the pvrolvsis of cellulos 
materials in a vacuum 
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5.2. Rates Calculated from Pyrolysis Data 


The data obtained in. the step pvrolvsis experi- 
nents at 280° (, shown in table 1, for cotton. hvdro- 
ellulose, and viscose ravon, were used in caleulating 
n terms of averages for 
ach individual step Because in the pyrolysis 
xyperimepts the degradation products were collected, 


veighed and analyzed, it was possible to calculate 


ates of volatilization for each fraction, as well as 


for the total volatiles These rates are shown plotted 


walnst percentage volatilization at 280° C in figures 
) 13, 14, and 15, for fractions V,,,. Ves, Vow, and 
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ty 
\ rt spectively : and in fivure 16 for the total vola- 
tiles The rates shown in figure 16 do not differ 
much from the Ay rates determined in the spring- 
balance apparatus and shown in figures 4, 5, and 6 
extrapolation of the rate curves for cotton and 


hvdrocellulose in figures 12, 13, 14, 15, and 16 leads 
to similar carbonization end points as were found in 
thr spring balance rate experiments In the case of 
viscose ravon, where pyrolysis at 280° ( 
only to 41 percent volatilization, the 
end poimt is not shown 


was carried 
carbonization 


6. Discussion 


the cotton used in this work had a lowe: 
and higher 
the hvadrocellulose 


Although 
ervstallinits degree of polymerization 
there is little difference in 
whether in regard to 
ratio ol products or to rates of degradation The 
crystalline forees seem to weaken at the pvrolvsis 
temperature, and, therefore, have little effeet on the 
results of p\ rolyvsis (As to the effeet of chain length, 
it was shown previously S| in the case of a number of 
polymers that long break up on heating 
during the initial stage of thermal degradation, even 
before volatilization occurs This tends to eliminate 
the effeet of chain length on pyrolysis, except where 


than 


ther behavior on pyrolysis, 


chaims 


polymers of very low degree of polymerization are 


involved Viscose ravon has low ervstallinitv and 


low degree of polymerization, and the type of its 
ervstalline structure is different from that of cotton 
and hvdrocellulose Results of pyrolysis of viscose 


ravon indicate a lower vield of tar, a higher vield of 
H.O, CO., and CO, and a faster rate of degradation 
than is the with cotton and hydrocellulose 
Further work will be required to determine the causes 
of the differences observed 

The presence of Na,CO, or NaCl in intimate con- 
tact 


Ciuse 


with the cellulose has a profound effect on the 


352 





results of pyrolysis of cellulose. The tat 


drastically reduced whereas the vield of HH. ('t) 
and CO is increased Although the aectua INI tig 
rates in the case of pure substances are low. rising 
to a maximum at about 15 to 20 percent volstilizg. 


tion, these rates im. the presence of Nat (). an 
NaCl start extremely fast and do not show any may. 
ma. The apparent initial rates are also much: high, 


in the presence of the salts than is the case w th thy 
untreated supstances 

As shown in table 7, the char vield is about ¢| 
sume for cotton and hvdrocellulose In the presen 
of salts, the char vield ts about doubled In the eas 
of viscose ravon the char vield is about 1.5 t 
much as in cotton or hvdrocellulose 

In the pyrolysis of untreated celluloses, up 
about 20 percent volatilization, the vield Ol tar rises 


ICS gs 


gvradually, whereas the reverse is true of the vields 9 
HoO. CO... and CO sevond about -0 percent Volatil. 
ization the vields are stabilized and reach constan 
values The data are insufficient to show whethy 


the same is true with regard to the salt-impre Triates 
celluloses 

\Mass-spectrometer, infrared, and e1 Voscople analy- 
sis of fraction V, all pomnt to the conelusion tha 
this fraction consists mainly of levoglucosan How. 
ever, there are definite Indteations that there 


Is pres- 
ent in fraction Vy... in addition to levoglucosa 
SOLE small amount of other compound Or cCOM- 
pounds TI ese Indications ure | i small] part ol 


fraction \, is insoluble in PLO, although all of it 
soluble in) methanol; (2) the infrared absorptio 
spectrum of fraction V,,, differs from that of 


clucosal In) having an absorption bani ato. , 


li VO- 


there are some differences in the mass spectra ol 
fraction \ and levoglucosan 

As shown tn table 7, the vield of fraction \ lal 
is 52 percent and 53 percent for cotton and hvdro- 


cellulose, respectively This means that iP pron 
mately | out of every 2 structural units in these cet- 
luloses vaporizes as a molecul ly 
viscose ravon, Which vields 30 percent tar, approxi 
mately | out of every 3 structural units appears 

the distillate as The next 
volatile fraction is HLO, which amounts to 21 percent 
for cotton and hvdrocellulose, and 29 percent fo 
After discounting the hvdrogen an 
oxvgen that vaporize from cellulose as parts ol levo- 
highest possible theoretical vield of 
other cellulose structural units woul 
from ana hvdrocellulose ani 
from ravon The difference he- 
tween the theoretical and actual vield of HO in th 
ps rolvsis of the celluloses Is 


levoglucosan 


levoglucosan larcest 


Viscose rayon 


vlucosan, the 
HO from the 
pereent 


he 27 ceotton 


3¢ percent Viscose 


aecounted for by th 


formation of CO, and CO, and by the presence ol 
Ho and O in the char That the chars obtained al 
the carbonization end points at 280)” ( 20% % lor 


and 33¢7 for 


contamed volatile 
whiel 


Viscose ravoth, as 
material Wis 


wert 


or cellulose 
table 7 
experiments in 


eotton 
vive in 


shown by these chars 








eld js 
(‘() 
Nitha 


rising 
tilize. 
ane 
Nay. 


1 }y, 


Varo- 
Prox - 


=“(' (Cl 





heate a quartz tube in a vacuum by means of a 
\Veker burner to about 1,000° C The loss in weight 
of the chars was roughly one-half for all the three 


eellu Os 


The results of pvrolvsis suggest the following 


poss! Thiet hanism lol the thermal degradation ol 
cellulose ln a vacuum \ dehvdration process takes 
place andomly along the chain \nothet process, 


h takes place Sl nultaneously, consists in thermal 


whic! 
sissions of the polymer chain resulting in the forma- 
rion of levoglucosan, H.O, CO., and CO. The C—O 
bonds in the chain are less stable thermally than the 
(CC bonds. This was shown by the present authors 
na ries of experiments on the pyrolysis ol poly- 
ethvl e,. and polvethvlen oxide The poly thvlene 
oxic vas found to have a lower thermal stability 
and n iower activation energy of thermal degradation 
than polvethvlene The various C—Q bonds in the 
chain are indicated in the structural formula (1) of 
the polymer as k, 1, m, and n 

Random SCISSIONS Of bonds k, I, orn res ilt ha more 
ol less complete breakdown of a part ol the chain, 


H.O, CO., CO, and C \ random scission of 


eld 


Hi | 
HoCon i 
| { {) (>) { a HH 
I] 
( ( ( 
oq] 1] 
HO ( 1 HO 
i] Ol 
P l)-y 1. DuA? 


the C—O m bond mas result in the formation of an 
oxygen linkage between carbons 1 and 6 by a process 
analogous to the formation. of levoelucosan from 
phenvl-8-p-glucoside MeCloskev and Coleman [9] 
have shown in formula (2) that this reaction takes 
place in two steps 

Whether in the pyrolysis of cellulose this 
formula (3), or in 
hvadrol 


process 


takes place In one step [10], two 


steps, is problematical. Pyrolysis is unlike 


vsis In an aqueous alkaline medium and a clos 
analogy between the two reactions is not justified 
The formation of levoglucosan units in’ cellulos 
pyrolysis takes place when scission of the C—O 
m bonds occur one. structural unit away from 
the ends formed through a random scission of a 
} () th) bond 

The presence of NacCO, and NaCl may have a 


catalytic effect on the dehydration and bond scissions 
taking place inside and outside the rings in the 
cellulose chain, resulting in increased vield of H,O 
CQO,, and CO at the expense of levoglucosan, and also 
In an increased rate of degradation of the celluloses 
particularly in the initial stages 
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Response Function of Thallium-Activated Sodium-lodide 
Scintillation Counters* 


Martin J. Berger and J. Doggett 


Measurements of gamma ravs with Nal 
related to the true ene rgv spectrum bv an integral equation Whose kernel 
probability that an incident photon of energy EF 
function has been caleulated by 
Results are 


is the 
The respol ~¢ 
energies from 0.279 to 4.45 Mey 


size from 0.25 (radius) by 0.5 inch Cength 


analysis of 50,000 photon histories sampled with the use 
Analvtieal corrections were 
and bremsstrahlung from the ervstal 


SEAC 


matic Computer 
radiation 


1. Introduction 


The response function of a Nal(Tl) ervstal spee- 
trometer relates the observed pulse height distribu- 
tion to the true gamma-ray spectrum. The primary 
determining this function is, and pre- 
sumably will remain, experimental. But a good 
theoretical understanding will contribute toward 
improving scintillation spectrometry. 

While the theory of the diffusion of gamma radia- 
tion in an infinite medium is now well developed, the 
ealculation of the response function is a boundary 
problem whose exact solution is bevond the scope of 
theory In its present state 

The subject of this article isa Monte Carlo caleula- 
tion of the response function, 1. e., an experiment 
carried out on paper. While this falls short of an 
analytical theory, it Is relatively simple to do, and 
has important advantages compared to a physical 
experiment. One is not hampered by the limited 
availability of monoenergetic sources nor by aecl- 
dental disturbing effects such as background radia- 
tion 

Exploratory Monte Carlo calculations of the 
response function have already been made by Camp- 
bell and Boyle ' for energies L=>6 Mev, and by Foote 
and Koch? for H24.45 Mev. An approximate 
analytical calculation for A=1.3 Mev was carried 
out by \laeder. \liiller, and Wintersteiger. These 
authors computed the effect of the first interaction 
of the incident photon in the cry stal exactly and made 
an elaborate estimate of the effect of multiple inter- 
This approximation limits the validity of 
their considerations to small crystals not much larger 
than one mean free path of the incident: radiation 


means of 


actions 





They confined their attention mostly to evlindrical 
crystals whose radius equals their length. In view 
*TI ‘ k Was supported t the I ted States Aton Ener (Com 
1.G. Campbelland A.J. F. Boyle, Australian J. Phys. 6, 171 (195 
RS. ind H. W. Koeh, Rev. Sei, Instr, 2, 746 (1954 
1 NI R. Miller, ar V. Winterst r, Hel Phy \ 27 4 
NI ] " n ely 
' ‘ ttw ’ i I 
Ty nt ‘ i \ } \ 
} fo ner na 


to 2.0 by 


Tl) eryvstals vield pulse height distributions 


response fTunetion 
will give rise to a pulse of size BE 
the Monte Carlo method for photons with 
presented for evlindrieal ervstals ranging in 
2.8 Y inches Thev are based on the 
of the Standards Eleetronie Auto- 
made for the escape of annihilation 


of the complexity of their “orders of scattering” 
calculation and of the approximations whose con- 
cannot be easily surveved, a «direct: re- 
calculation of their results seemed of value. More- 
over, one would also like to know more about the 
influence of the crystal shape on the 
function. 

Monte Carlo calculations have been made of the 
response functions of evlindrical crystals of diverse 
shapes ranging in size from small to very large for 
radiation incident with energies of 0.279, 0.661, 1.17 
1.33, 2.62, and 4.45 Mev 

The absorption ana ¢ ‘ompton scattering of photons 
in the erystals were calculated by random sampling, 
the computations being carried out on the SEAC 
Approximate analytical corrections were applied to 
the Monte Carlo results to account for the escape of 
secondary bremsstrahlung and annihilation radiation 
from the erystals.° These corrections were quite 
small at 1.17 and 1.33 Mev, but appreciable at 2.62 
and 4.45 Mev. Because of the approximate nature 
of the corrections, the most important and accurate 
parts of this work are those pertaining to the four 
lowest energies where the corrections were negligible 
Most of the results are for collimated sources, while 
some sample calculations pertain to broad beams, 
point-isotropic sources, and off-axis collimated sources, 

The main results of this investigation are: (1) The 
confirmation and expansion of the results of Maeder 
et al. for small ervstals; (2) the extension to large 
erystals in which multiple interactions are impor- 
tant and to higher energies; 3) last but not least, 
the creation of a SEAC code by means of which re- 
sponse functions for source and crystal geometries 
peculiar to a given experimental situation can be 
readily computed with a modicum of effort 


sequences 


response 


for this pproximate procedure was the ten porary not 
' 


‘ ison 
labilitv of SEAC during the latter part of this investication It would 





have heen he rt rry through the entire calculation on the automatic con 
puter rhis is particularly sy for quanta resulting from the annihilatior 
positror oduced in first interaction of the primary photon The source 
but he phot within the ervstal can easily be levilated, and 
ut yuent h ryt t t ! of the exist SEA cnele We ( 
nge the SEA “le i he that it can handle the escape of 
n \ lintior 1 t ‘ } e44 Nie In the 
t t ‘ ( er t ese esu ind 4.4 
Ne ce ene el bee ‘ rte t is] 








2. Method of Calculation 
1 Formulation of the Problem 


ponse function can be expressed in the form 


Lik [: is the 

nitial energy / will deliver to the ervstal 

sult of one or more interactions—an amount of en- 

between A’ and k’+dk’: and q(kh’ kdb” is 

e probal that upon delivery of energy /’ to 

Stal, the light output as amplified by the 

photomultiplier will indicate an apparent 

heiht between ’’ and hk’ lat 
experimental fact that 


probability that a photon 


us the 


? \ 
yrgity\ 


wrey (pulse | 


2rak’ eXp hk’ —E"’ ak’ 2 
he value of the parameter @ depends on the physical 
characteristics of the ervstal-photomultiplier system 
We shall consider this parameter as determined ex- 
perimentally, and concentrate on the caleulation of 
the loss function L( ik’ 

To clarify the physical pieture, it is useful to break 
the loss function up into a number of component 


parts, writing 


L(E,E)=Y (EF) K(E, E’)+ p(£)s( EE’ 3 


(/) is the efficiency, 1. e., the probability that a 
photon incident with energy Ie will have at least one 
in the ervstal 


Interaction (seattering or absorption 


For i collimated source 
Y()=1—exp [—n(E)L 


where L is the length of the cry stal. and ul ft the 
narrow beam linear absorption coefficient 

A( hk dk’ is the probability that an interacting 
from the erystal after one ot 
hore scatterings, leaving behind an 
gmount of energy between /’ and h’+dk The 
funetion p(/e) represents the probability. that) an 
interacting photon will be absorbed (possibly after 
a number of Compton scatterings), while the delta 
funetion 6(4—") indicates that the entire energy 
fy is delivered to the erystal in a photon history 
in absorption. We note the normatliza- 


photon \ il] emerge 
Compton 


terminating 


Lion ‘ 
dE’ K(E, EF’) plik 5 


For a wide range of conditions the shape of the loss 
funetion depends largely on the ratio of the numbers 
of absorbed aad emergent photons, and only rather 
K(k Ee’ In a rough approxima- 


Insensitively on 





characterized py 


shall « L the 


Lion, the shape ean therefore be 


following Maeder— we 
In experimental terms, 


p(k), which 


photofract 1On 


area under the “photopeak”’ ol the it] se 
height distribution 

WE: , 

urea unde! the entire pulse heieht dist! tiny 

' 

Knowledge of the photofraction is useful | Othe 

cConnectlons Thus if one wants to use a cry al 10) 


counting photons in. the 
eround, il That \ be desirable to count only t| lapoy 
pulses resulting Prom complete photon abso Dtioy 


The required “photoefliciency” of the cry tal is 
) , 


hy) p(k Counting of photons in the ‘photo. 
peak” only is particularly useful for disentangline 
complex spectra containing several lines 
2.2. Random Sampling 
Key the calculation of the loss funetion, i direc 


unalog method 


stochastic paralleling the plivsica 
processes i all respects, has been chosen Seu tilla- 
tion ervstals are verv efficient detectors; hence thy 
efheiencs ol an analog Nlonte Carlo calculation S 
also high, and the computational cost of statistics 
and analytical refinements would be out of proportioy 
to the possible gain in efficienes 

The caleulation proceeds along the following lines 
We start a photon al 
ervstal surface with specified energy and direction 
The position of the first: interaction 
ys if the photon moved im an infinite scattering 


il specified Position on thi 


is determined 
medium: the distance traveled to the point of inter 
netion is distributed exponentially, the mean of the 
distribution being the inverse of the total 
beam absorption coefficient If the lniter- 


Hhear 
Harrow 
action occurred outside the region occupied by th 
crystal, the photon is recorded as having had 

interaction, and another photon history is beou 
If the interaction ts inside, we determine by randot 
sampling whether it is an absorption or Comptor 
scattering, the relative probabilities for these contin- 
CeHCLES beme proportional to the respective Cross 


sections ln ease of an absorption, the photo 
history is terminated, and another begun lin ease 
the new energy and direction of th 
photon after the collision are sampled from th 
Klein-Nishina distribution, whereupon the loeati 
of the next interaction is sampled in the same mann 
as before. If this 
ervstal, the photon is considered to have es auped 
and the ene rev left behind in the ervstal is recorded 


proceed iis befor 


of a seattering, 


interaction Occurs outside the 


If the Interaction was inside, we 
deciding whether it was an absorption or scattering 
ete The procedure Is repeated until thre photo 
has either left the ervstal, ol has reached an energ' 


RK 95 i 


W. EK I 
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low than 50 kev, in which case the photon is 
The introduced by 
this toff is negletble be ause of the steep rise of the 


Nal tt 


eons lered as absorbed bias 

pho lecetric absorption Cross section for 

erie es below 50 kev 

T detailed equations for carrying out the 

ean ne procedure have been presented elsewhere. 

together with an account of the numerical manipula- 
eq ured to adapt the 


Various Cross sections to 


on Itaution ona hich speed compute 
The required absorption coefficients were obtained 
7 hie tabulation of G. White. 


2.3. Machine Computation 


he SAC program the caleulatior 


fo) e simultaneous consideration of nine erystals 


Was set up 
of different sizes Kor the sike of convenience these 
sizes were chosen so that the ervstals formed a nested 
each ervstal containing completely the next 
Whena photon in the random sampling 


sequ ce 
small Prone 
xperiment eseaped from a erystal, the historv was 
continued for all crystals of larger size. This pro- 
dure resulted in great economy of computing time 
Bevo d that, it the statistical efficiency ol 
comparisons ol the response lunetions of ervstals of 
different shapes 


The output ol 


hnereases 


a SEAC ecale 
Information for each of nine crystals 

a. The efficiency (4 This information is ina 
edundant Ssinee it can also be easily obtained 
vou direct analytical calculation eq } sut a 
omparison analytical Nionte Carlo 
ferences pros ides a usel il che kk thrat the compute 


tution consists of the 


followine 
between nnd 
operation Was error-Lree 


b. The spectrum (in histogram form 
packets supplied to thre ervstal hyy 


of the enere’y 


photons whi hy 


I: K(k f 


ce. The photofraction }) I: 

lt was conventent to generate and process photon 
stories In groups ot 1.000 Kor photon S with an 
do to 40 min 


energies the running time 


itinl enerev of GEO key 


of SEAC time: for 
ws a little longer 


this took about 
hich ! 


2.4. Escape of Secondary Radiation 


Nlonte Carlo caleulations it was assumed 


entire enere’y 


iT the 
that the released im the 


scatterings or final absorption of the incident photons 


SUCCESSIVE 


he Bs 56 ! / R P24 


remained 


in the crystal Let rally this is the case 
only ni low energies At enereies above the threshold 
for pair produ tion, annihilation quanta may escape 
from the crystal: additional energy Cin he lost I) the 
form of bremsstrahlung. In the energy region of 
up to 4.45 Mev) the effect of bremsstrahlung 
was found to be minor compared to that of annihila 
tion radiation Not only is) the lattes 
quantitatively, but it also affeets the shape of the 
function in a more distinet fashion YIVIng 
rise to two additional pe aks corresponding Lo thie 


Interest 
ereatel 
respolise 


2 annihilation 


. ibsequent to the absorption Ol a primary 


complete escape ol l o1 quanta 
photon 

The appropriate corrections were made according 
to the approximate outlined in’ the 
appendix. The photofraction p(/’) was 


both for annihilation radiation and bremsstrahlung 
In ord rto avoid CX\COSSLV COMPULAtiOn, the detailed 


proce lure 
eorrect a 


shape of the response function was corrected only for 
annihilation radiation 
effects of bremsstrahlung was made bv a subsequent 
renormalization of the funetion AChE’ 


and an adjustment for the 


2.5. Gaussian Broadening of the Response Function 


| 2), (5 and (4) that the 


It follows from eq 
response Tunction 


Y » 


\ Irak 


REE” 

| si 
— @Xp ID 
fe’ Pale’ dk’ > 7 


exp {—(E”—E)! 20k: 





H( BE E’)+GQkh bk” 


For numerical applications the Was 


0.00146 Fe(\Mev corresponding to a 


half-maximum W=0.09,/. This is a 
typical value in close agreement with experimental 
data of Foote and Koch (see footnote 2 

In the evaluation of /7(/2.4°") one must take into 
account that AC, /’) is available only In histogram 
form 


parameter a 


ehosen to he 


width at 
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mead that the k's have been corrected if 
for escape of secondary radiation. We 
hem to form a row matrix &(/ koh 


We similarly replace //(/,/°’) by a row matrix 


/ ho) } where h, (ke Mi hR dE" 


convolution of ACA’) with a Gaussian ts 
replaced DV a matrix multiplication 

hi Tk. S 

where the matrix element 7 a kind of transition 


probabilits Is vivel by 








3. Results 


3.1. Efficiency and Photofraction for 
Radiation 


Colli: ated 


(rvstals of such sizes as are either readily 
able commercially, or are in use in this laboratory. 
were selected for computation In table 1, the « 
efficiencies are listed for collimated radiation of 
various energies incident along the axis of thy 
drical crystals, computed according to eq (4 Che 
figures in parentheses indicate the values that wer 
obtained by the Monte Carlo caleulations. The 
agreement is excellent 

In table 2, the photofraction p( 
is the corrected value taking into account the ¢ scape 
of secondary radiation. The amount Ap which had 
to be subtracted from the raw Monte Carlo result 
Po in order to make the correction ts indicated jn 

The indicated 
computed as follows 


0 B BP 10) 
VN 


where .V is the total number of photon histories con- 


is given. This 


errors are. standard 


The 


parentheses 
deviations 
deviation of py is 


standard 


tributing to the determination of p Nis also 

shown in tables 1 and 2 The correction Ap was 

} t ) . (' ‘VE assumed to have a standard deviation (1/10)4) 
2 _ 2 ; This is a nominal figure representing an “educated 29 

CuUeSS 4 The total standard deviation is 
The exact location of a pulse within an interval of | 
size & 32 is not known from the Monte Carlo eal- T= |% 1/10Ap 0 
culation We assume a uniform distribution, which 
is the reason for the average with respect: to 7 in M ; S| 
| } I ( ff ; ‘) 
g Ni ( 
\ ‘ ( 
kk a? 
\l 
' \ “Ih 
s 
ery 
hs el 
~ ‘ rie 
| 
nA=0 











CO])- 


also 
Y Was 


‘ated 


pends on L 


TABLE 2 


3.2. An Empirical Formula for the Photofraction 


It has already bee 
for small ervstals at 


nh noted \Iaeder et al . that 


of O51] \ley the 

funetion of y LR? where L 
ervstal length and AR the radius The results 
that this a true for large 
In the entire from 0.279 


hy 
ahh Cle rg\ 
phototraction 


s thy 


Is a SMOOT! 


+) 
Ol table 2 si 





ow st) holds 
rvstals and 


to 4.45 Mev 


\loreovel if 


energy rang 
he cle pence nee of the photofraction on 
erystal size can be accurately represented bya 
simple empirical formula, under conditions where the 
tion due to the ese upe of secondary radiation 
bsent small (i. e., at 0.279, 0.661, and 1.17 
\lev for all thr and at 1.33 for all but the 
smallest ervstal considered 

The formula ts 


} 
ervstais 


here L and R are the length and radius of the 
ervstal With the parameters A, B, and © listed 
1 table 3, formula (12) represents the values of 
(/-) given in table 2 with an accuracy that is 1 to 2 
percent in most cases, occasionally 3 to 4 percent, 
nd in the worst cases 5 to 7 percent These devia- 
ons are presumably due to the facet that p fe) de- 


and R 
rough \ LR 


in a more complicated way than 





Photo ym pit 
{ | 
\ 
PABLI / } p ) 
; } “i ; 7 
I \(/ 3(/ C'(} 
\/ 

0. 279 0. 47 0. 514 QO. 65S 
Hol O40 S4] 21% 
17 Xt) yd 124 
| IS5 Y42 120 


Formula 
albedo for a 


and R 


number 


When L 


12) also predicts the correct 
semi-infinite medium 


Lh l number albedo 13 


When we extrapolate in’ the other direction 


L QOand R-0 
p(Ok 


Lh 14 


where up, uw is the ratio of the photoelectric to the 
total attenuation coefficient Thus, the three 
parameters A, B, and C, two are fixed by (13) and 
14), and only ©, a quantity with the dimensions 
of an attenuation coefficient, ts adjustable This 


makes it rather remarkable that a simple formula 
like 12 


of 


Wol ks 
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3.3. Shape of the Response Fu 
Radiation 


The response functions for all the ervstals listed 


in tables 1 and 2 except No. 2? 


la and 1.b for ineident energies of 0.279. 0.661. 


1.17, 1.33, 2.62, and 4.45 Me 


response function exclusive of the photopeak, i. e energy being covered by 32 intervals of equa 
~4 
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photopeaks are indicated by arrows 


The Gaussian broadening has been cal 
according to the preseription given in sectio 
It should be mentioned that the shape of the res 
function is rather insensitive to the magnitude of , 
\ For the sake of | The results are presented in histogram forn 
clarity ol presentation we have plotted only the ; 


ure shown inh figures 
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‘| histograms are normalized absolutely so that | contribution to the scintillation spectrum due to thi 
thie correspond to one photon incident on a escape of secondary radiation In other words. they 
ervsial. The area under each histogram is thus | represent pulses which would have ended up in the 
Vik [1—p(-)| In order to supply the missing | photopeak if the absorptions had not been spurious 
photopeak, one need merely add a Gaussian with | It can be seen that at sufficiently high souree energies 
thie lesired standard deviation ef eg 7 and the response functions have two peaks nt energies 
mul plied by | welght factor ) Ik p Ie ~ a, kk nic and 2 me below the incident eneryy ke. which 

and Ie) ean be found in tables 1 and 2 ire associated with the escape of 1 or 2 annihilation 
The shaded parts of the histograms indicate the | quanta 
It can be seen that the shapes Ol the respons 
functions for various crystal sizes, and even for 
different ehergies, are on the whole rather similat 
of 
NO.5 b No.4 , »& | nos , § © NOL! yp § 
R= 0.75" f R= 0.875" H R= 2.5 H R=2.5 1h 
L«2.§ ni = 2.0" iy i 4.0" L= 9,0" i] 
he} mi) “* #}) 8. iI 
wee abs HH 
rf ve A ee ir 
be id Hi UY r 6 
3 E= 445 Mey tr } | ve p 
10 aes nt + {| i U a 4 
~ $4 Se ie oa” Hi M } 
7 A pe pe Sine onl i n 
Sl pu L — Jha f] ru _ i elr g 
o Ih 4 - hh H LJ Ze J ] 
nm in] " irs . Ay r -) ot ‘] 
HL, f L ¢@ } mr | 
2 E:262Mev pg & ml me | A g 0 
a-4 raf ia u a 
i aa en dl 7 rd 
= — a, ed | - 
— lod w S| “Ls fs) 
} W Uj 
u q if T 
E*1.33 Mev 
rm L 
<a —_ aw) rv r 
4 are al ~ aw I di iF 
io + ™ LF L \ 5 ol |} # -_ tH 
LF | Y - a t — . Hi 
h. L 7 —_ 
E= 1.17 Mev - - 
a=2 we, ae oda 
ee 4 — 4 * “al ‘7 ie Ar " 
O .™~ | “Ly “es Le } *: pu L - ~J L 
L i! ne er om 
J — L Lu ; Ll) 
= E+ 0.66! Mev L ! 
rt a=| SL ‘7 rr) 
™ 4 a i = . 7 a) - 4 
in bite cats = ee | L 
Ml P) 7 + ~ 
oO } } eine al ‘bn, = \ > n ++ 
s) 7 7 4 5 | 
Oo 7 ane m J 5 
_ Pr &, 
Ss L } ; 7 
<2  €*0.279 Me L 
z a-0O 7 ale 
—, ” a u 
rm 2 L uy ~ + 7 
Le 6G ia \ m 
Sy — 
3 t 5 uW 
O t | } ; 
L L 
3 L . 
4 7 7 j ] 4 ] 
10 ‘ LU J 1 jaewere | J — 41 
0 8 6 24 32 0 8 16 24 32 0 8 16 24 32 0 8 16 24 © 
32 é 
Pe Pint ; Nal | } if } ng tt ( '¢ 





361 





( 


iil 


al fluetuations§ of 


is ey idenced 
are much 
energies It 


ls 


al irregularities 


by 


the 
the 
worse for high 


Mionte Carlo 
irregularities of 
than for 

that quite 


is noticeable 


occur 


for 


histograms 





taining to the same energy but different ervstals 
Sphenomenon is due to the use ol the S2aine Froups 
all the ervstals 


) photon histories for 


3.4. Other Source Geometries 


Monte Carlo calculations of the response function 


lol ny desired 


easily be 


source geometry can 
necomptlished tl rough the proper choiee of the initial 
eonditions, photon positions ana clirections, in the 


sampling Sample results of such 


O66] 


random process 


caleulations for Nev radiation are presented 


in figures 2, 3, and 4 
a. Broad-Beam Radiation 
n figure 2 the ratio Pua pis plotted against 7 PR, 
Where Ppp bs thr photofraction for a broad ( viindrical 
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pis the 
fraction for a collimated source, and FP is the « 


beam with cross-sectional radius ;. 


radius 


large: No. 6, medium-sized: and No. 9, small 


b. Point-Isotropic Source 


In figure 3, the ratio prs p is plotted ag 
tan”! RD where pys is the photofraction for a 

Isotropic located at a distance J) fron 
crystal on the evlinder axis. The curve perta 
crystal No. 6. The end-point of this curve 
tan 'R D0 was obtained from the COrrespol 


broad-beam result 


source 


c. Off-Axis Collimated Beam 


the shape of the response funet 
pencils of radiation ini 
from the crystal axis 


No. ] 


In figure 4, 
shown for collimated 
at various distances s 
crystals No. 4 and 


Three curves are shown, for ervstals \ 


OO. 
Stal 


tinst 
omt- 
the 


is t 


for 


dling 


(1) Is 
dent 
fo 


The normalization and 


other details are the same as in figure 1: the photo- 


fractions are also shown 
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OLo- 4. Discussion ean equally well be applied to our results. It indi- 
Stal cates fair agreement between theory and experiment, 
VO. 9, l. Comparison With Other Calculations but there is a consistent tendency for the experi- 
Our results for the value of the photofraction for | mental photofractions to be somewhat lower than 
ollimated radiation incident on small crystals are | the theoretical values. This diserepaney may be 
n cood agreement with those obtained by Maeder reduced in part if one takes into account energy 

tinst Pet al For example, at an energy of 0.661 Mey losses due to electron escape from the ery stal 
oint. Jone finds: Measurements of Foote and Koch. see footnote 2 
the indicate for collimated radiation incident on a large 
ns te oe oe ervstal (?=2.5 in., L $ in.) photofractions p=0.8 
San at 0.661 Mev, and 0.64 for cobalt-60 radiation, while 
ndine : the corresponding theoretical values are p= 0.821, and 

e 8 iD 5 } 0.675 respectively." 

Further results for a large erystal (/?=2 = in., 
\iaeder 0.246 O.368 Oc. 418 O. 471 0. 507) OL 546 L fin have recently been reported by Kreger See 
! 2 245 Sia $42 1S] 50S o42 footnote 6). His incident radiation was a broad 
On is beam with cross-sectional radius r=0.25 Rh. Hence 
dent the comparison must be made for the broad-beam 
fo At 0.279 Mev and 1.17 Mev the agreement is photofraction Pas. We have so far only calculated 


and Jalmost as good (with the exception of a 15 percent | this correction at an energy of 0.661 Mev, but esti- 
hoto- iscrepaney for the smallest crystal No. 9) at 
1.17 Mev 

No other results comparable to our photofractions 
for large ervstals No. 1, 2, and 3 exist. Maeder et 


mate that with an error of only a few percent the 
same fractional correction can also be applied at the 
other energies of Kreger’s experiment. 


* Ta. did, however, attempt to compute the photo- — 
fraction for the limiting case of an infinitely large 
rvstal (i. e., a semi-infinite medium), calculating - 
the first two interactions of the incident photon 0.209 0. 661 ly 2 2. OR 
cactly. and estimating the effect of the subsequent 
nteractions. We have calculated the photofraction » (iineen 0.93 0.725 0.57 O34 0.47 
for this case by extrapolation based on the empirical p (present res ilts ) 74 oOo aS 0) 
ormula (12) and arrive at the following: comparison: 
Knerg Vey The agreement is seen to be quite good at 0.661 Mev, 
mn at where the exact value of the correction is known, 
4 > . 7 and fairly good elsewhere 
|| 0.265) 0. O61 Ld Summing he entire experimental evidence, 
L uimnming up th nti xpert 
both for small and large crystals, we can state the 
= p’ Monte Carlo O. O75 0. O40 0. USO experimental photofraction is invariably helow the 
Maedet Wo Sib S02 theoretical values, by a smaller or greater amount. 
~ One is led to the conelusion that in the experiments 
54 there is always background radiation present (such 
L || The “orders of scattering” approach, in the limit | as back-seattered radiation from the radiation source 
f very large crystals, agrees with the more exact | and from the material surrounding the detector 
aleulation only at 0.279 Mev, where multiple inter- | whieh would be counted in the tail end of the pulse- 
74 rtions are unlikely because of the large amma-rayv height spectrum and would thus tend to depress the 
bsorption cross section. but leads to a sienifieant value of the photofraction. 
nderestimate of the photofraction ut higher energies 
= Qur albedo values | p are in vood agreement with 
] w results of a calculation by Hayward and Hubble The authors thank Dr. Evans Hayward for several 
Vote that the albedo is greater at 0.661 Mev than at | stimulating and enlightening discussions, and Miss 
279 or 1.17 Mev. There is also experimental evi- | Mary Orr for her assistance with the hand com- 
1] ence for the oecurrence of this maximum at an | putations. 
L termediate energy 
ot § The pictures of the response funetion given by 5. Appendix: Escape of Secondary 
. Maeder et al. are similar in appearanee tO Ours, but Radiation 
ot detailed enough for an accurate comparison Ih) 
gard to shape 5.1. Effect on the Photofraction 
4.2. Comparison With Experimental Results To estimate the effect of the escape of secondary 


radiation on the value of the photofraction p, we must 
: ! > PL y . *}° . ‘ 
Fon smal] crystals HT detailed comparison with CA- determine the probability that an absorption Is 
erinents has been made by Maeder et al., which 

, K. Lidén and N, Starfelt, Arkiv Fysik 7, 428 (1954) : 

EH iJ. Hubb a) R 93.4 4 +40 radiation. ae j +) 


} nhotofractior } } 
WwW \ | } . ' ne . , ~ | 
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e., that it results in, or is preceded by, 
ap of at least secondary photon from 
tal (with an enerey greater than 50 kev)." 

We shall ig the displacement of charged 
the ervystal." Bremsstrahlung and 
hilation quanta are thus assumed to be emitted 
he pomtl of absorption of the primary photon 
\loreovet the emission to be 
| I ut the 
ytions which occur in the first interaction of the 


one 


nore 


is assumed isotropic 


tention is econeentrated = on spurious 
prior to any (‘ompton seutter- 
unimportant 


other imter- 


nil photor 1. @ 
nd return later to the relatively 
from 


nee 

rs il 
resulting 
™~ condary radiation is therefore considered 
from an exponentially distributed line 
on the assumption 


spurious ubsorptions 
actions 
ns originating 
the axis of the ervstal 

radiation is collimated 
notation will be used for the various 
photoelectric absorp- 


SOs’ dt) 
that thre incident 
The following 


coefhieients 


attenuation UpH 

tion coefficient, up =pair production absorption co- 

efficent, we =-Compton seattering attenuation coeffi- 

cient J UPH Up Lic total attenuation coefficient 
AD hilation Radiation If an absorption of a 


primary photon of energy / results in pair formation 
and subsequent annihilation of the positron, the 
probability that at least one of the two annihilation 


quanta will eseape from a evlindrieal erystal (with 


radius 7? and length Z) is given, approximately, by 
the following expression 

A(T [ ) , 
) ; (id¢ 
f eXp lu ip / . 

d cos @. ] | ( ] ( Pip? LS 
where me and and s are the distances 


traveled in the ervstal by the two annihilation quanta 
emitted in Opposite directions 





_, 
If 0<r<L2 
Cos COs « S rcosé N IF r) cosé 
eos COS 1 OO COsS 4 S rcose N»9 Pe /sin ‘ul 
“cos §< cos 4 S sin @ No PR sin 6 
If Li2<r<b 
COs 6, cos @< | S rcosé N. If r) cos @ 
cos 7 cos # \cos 6, Ss I sin 6 S L r) cos@ 
(“cos 6 < cos 6 N PR sin @ . P sin é 
4 
16 
where 
R i R 
Cos 7 | ( ) » COS 7 | ( ) 
L—s : J 
| ( 
LOM electrons in Nal are approx ely 0.02, 
wi Hence only at the highest ener $45 Mey 


two anni ation 


pe and p¥ are photofractions for the 


quanta, which—we believe—can be accurately agp. 
proximated by the photofraction p for an energy 9 
0.511 Mev, as obtained by interpolation from (able9 


Bremsstrahlung The probability of esc) of 


bremsstrahlung quantum of energy /’ can be ey. 
pressed by a formula similar to (15 
Pik’ wf , dare 

| CXp ul fe yh e 

, 

: dcosé - ¢ ( | ) fk 
where we ugaihn set y* (fk p hk and tl at 
length S Is defined iis follows 

I 5 
cos 6 cos 8< 1: S 
: cos 
R 
Cos cos 4 cOs F.° \ > g 
Sind : 
r—L 
COs COs f° S 
COS OJ 





cos 6 and cos & have the same meaning as abov 
Let 4) Kk, ke’ dk’ denote the probability, 
path length, for an electron or positron of ene rev | 
to produce FT| photon in) the eChergy interval k 
k’+dk’), 


per unit path length of an electron or positron du 


to 1lonrIZINg collisions The probability of ahh escap 
of a bremsstrahlung quantum produced by 
electron or positron of energy Ie in the co Ise ¢ 


slow ne dow Ih is 


PAR 


where R(/-) is the electron o1 positron range, al 
E 50 kev is the cutoff 

Let Wihjdk be the probability that 
members of the pair has acquired initially a_ fra 
tion between & and |} dk of the total availabl 


kinetic energy The probability of the escape | 


+} 


one ol 


at least one photon in the form of bremsstahtung 


emitted by a pair of total energy fe while slown 


dow no is 


Pik dkWk)\|PARKE)+P ARK 


PARE AP AE kh 2 


Hl. Bethe W. Heith I K : \| 146. 5 
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per unit 


; and let dk ds be the average energy loss 
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ator 
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ible? 
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slow I 





The above formula assumes that each member of 
he ir cah emit at most one photon with energy 
gren than 50 kev in the course of slowing down. 
We eve that the effects of plural photon produc- 
tol in be safely neglected in the energy range of 
nite st here 

Combined Effect. The probability that the ab- 
sorption of a primary photon of energy / in a first 


combined 
and annihilation radiation 


nternetion will be spurious due to the 


effects of bremsstrahlung 


In evaluating this formula we found that the nee- 
ect of bremsstrahlung (i. e., setting P?. and P,—0 
would lower Py by only a few percent (a maximum 
of ~ 11°) for ervstal No. 9 at 4.45 Mey The prob- 
abilitv, per incident photon, of a spurious absorption 
n the first interaction ts 


PAZ VE) w(E)| | epn (1) + up (Ee 


Vext we must determine a corresponding probability 
of spurious absorpcion in a second interaction, One 
ontribution to this probability is similar in origin to 
21), and can be written 


- - 
PAE Vik) ult Zz = F hse f 
e } LL Ik 
KEE QE PRE’), 233 
vhere ki FF’) is the Kleim-Nishina differential co- 


ficient for Compton seattering with energy change 
from i to Ek’: (J is the probabilits that the seattered 
ncident photon will not eseape from the ervstal with- 
uit further interaction *'; and Ps fe) is the probabil- 
tv that a given absorption event is spurious, and is 
inalogous to af Ie but should take into account 
the facet that the source of secondary radiation is no 
onger concentrated on the ervstal axis. We have 
gnored this circumstance, because the evaluation of 
in exact expression for P¥(/2) would be exceeding] 
omplicated and have set re P, Ie It is 
oped that this will introduce only a minor error, 
or thi For small ervstals, such a 
zeometrical error might be serious, but a second inter- 
etion is quite rare; for large ervstals, in which multi- 
common, the geometrical 


following reasons 


We interactions are more 


anv ense, the 
phototraction 


error is much less significant In 
contribution of 2? to the change in the 
will be small, so that a rough-and-ready treatment ts 
Ih ordet 

There is also the possibility that 
taking place in a second interaction is spurious be 
cnuuse it Was preceded by the brems 
strahlung quantum from a Compton recoil electron 
The probability for this to happen is 


an absorption 


escape of a 


.7 - 
Pik Vik) ulk dE’ = ~ rf 
Ji u(t: 
KEE QE )P (EE 24 
Thus a spurious interaction will oceur with 
probabilits 
P(E)=P,(E)+PH(E)—P(E)PHE 25 


per incident photon in a second interaction; and with 
probability 


Pi(b)=P(E)+P(E)—PS(E)PACE 26) 


in a first or second interaction \t the highest 
energy and for the largest ervstal considered, the 
value of Py was only 15 percent of P;, and generally 
it was much smaller. Subsequent interactions will 
lower the photon energy further and make spurious 
absorptions even more unlikely. Hence we felt jus- 
tified in ending our calculations with the effects of 
Interaction The change in the photo- 


fraction is thus 


the second 


Ap PAE)Y(E). (27 
5.2. Effect on the Shape of the Response Function 


In considering the detailed shape of the response 
function we have calculated only the effects of an- 
nihilation radiation from a first interaction. Other 
escaping secondary radiation was taken into account 
only insofar as it shifted pulses from the photopeak 
into tail of the response funetion (i. e., by an appro- 
priate renormalization of A(/, 2” 

None, or one, or both of the annihilation quanta 
resulting from a pair production can escape from the 
crystal without further interaction. The probabil- 
ities for obtained by 
replacing the integrand in the inner integral in (15 
by the expressions 


these contingencies can be 


tl ho escape / ¢ "i" é Y 
ly one escape e- ee ( 
2S) 
‘ 
¢ Iwo escapes ( a a J 
> 
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y hk K me, 


e terms; 


two 


contribution 


2me 


b), #’ is the 


ee ording 


to 


fe’ 


Is 


sum of 


aT) 


left im the crystal is the 
of them are each distributed 


p mc 


constant 


the 


two 


other 


6 hk ‘ 


al 


Trig IQ 


mount of 


terms, one dis- 


constant 
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equal to k—2me*. In ease (c) k’=kH—2me The 
correction of the response function was performed 
by shifting the appropriate number of pulses froy 
the photopeak to the ‘tail’ of the pulse leigh 
spectrum, distributing them according to t} (Is. 
tribution functions described above weighted 
cording to (28 


iat’. 
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